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ABSTRACT 



We present K band spectra for 133 nearby (d < 33 parsecs) M dwarfs, in- 
cluding 18 M dwarfs with reliable metallicity estimates (as inferred from an FGK 
type companion), 11 M dwarf planet hosts, more than 2/3 of the M dwarfs in 
the Northern 8 pc sample, and several M dwarfs from the LSPM catalog. From 
these spectra, we measure equivalent widths of the Ca and Na lines, and a spec- 
tral index quantifying the absorption due to H 2 opacity (the H 2 0-K2 index). 
Using empirical spectral types standards and synthetic models, we calibrate the 
H 2 0-K2 index as an indicator of an M dwarf's spectral type and effective tem- 
perature. We also present a revised relationship that estimates the [Fe/H] and 
[M/H] metallicities of M dwarfs from their Na I, Ca I, and H 2 0-K2 measure- 
ments. Comparisons to model atmosphere provide a qualitative validation of our 
approach, but also reveal an overall offset between the atomic line strengths pre- 
dicted by models as compared to actual observations. Our metallicity estimates 
also reproduce expected correlations with Galactic space motions and Ha emis- 
sion line strengths, and return statistically identical metallicities for M dwarfs 
within a common multiple system. Finally, we find systematic residuals between 
our H 2 0-based spectral types and those derived from optical spectral features 
with previously known sensitivity to stellar metallicity, such as TiO, and identify 
the CaHl index as a promising optical index for diagnosing the metallicities of 
near-solar M dwarfs. 

Subject headings: stars: late-type, stars: fundamental parameters, stars: abundances 
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Introduction 



Contrary to their dimness, M dwarf stars hold significant promise for illuminating the 
processes that govern the formation and evolution of stars, planets, and t he Milky Way 



M dwarfs possess masses of 0.6 M > M* > 0.08 M (jPelfosse et al 



20001), straddling the 



peak of the stellar i nitial mass function and dominating stellar populations by number 



( iBastian et al. 



20101 ) . With main sequence lifetimes that exceed a Hubble time, Galactic 



M dwarfs also represent a complete arch eological record of the chemical evolution and star 



formation history of the Milky Way (e.g. 



Bochanski et al. 



20071 ). Moreover, M dwarfs are of 



great interest as potential exoplanet host stars, as all of a planet's observable signals will be 



significantly easier to de tect if it orbits an M dwar : 
around a G dwarf (e.g., iNutzman fc Charbonneau 



com pared to a similar planet in orbit 



20081 ). 



Calibrating the fundamental parameters of M dwarfs, however, is a difficult challenge 
from both an observational and theoretical perspective. Empirical measurements of M dwarf 
masses, luminosities, temperatures, and radii are anchored primarily by the information 
extracted from the orbits o f M dwarf binaries , either the very rare eclipsing systems 



>15 systems known; e.g., iKraus et al.l 1201 ll ; llrwin et al. 



201 lh or 



somewhat less in: 



Martinache et al. 



'orma tive) spectroscopic/astrometric systems (e.g., 



ess rare (but a 


so 


Shkolnik et al. 


2010 



20071 ) . Overcoming the intrinsic faintness of these low- luminosity systems 
however, makes their characterization a taxing observational challenge, and analyses of 
known binaries reveal systematic offsets in inferred temp eratures and radii that correlate 



with both orbital period and tracers of magnetic activity (IRibas 



2006 



Lopez-Morales 



20071), 



suggesting that m ore wide binaries a re needed to accurately infer the parameters typical 



isolated M dwarfs ( IKraus et al. 



20111 ) . Theoretical constraints on M dwarf parameters have 
been similarly difficult to achieve: accurate l y modeling the deep convective zones in M 



dwarf interiors ( jMullan fc MacDonald 



2001 



Browning 



20081 ) and the mixture of molecules 
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(Tsuji et al. 


1996; 


Allard et al. 


2000) 



requires significant computational resources, as well as extensive databases of oscillator 
strengths and opacities measured from laboratory experiments. Confronting theoretical 
predictions with empirical measurements have also identified significant offsets even for 
isolated field stars: empirical effective temperature scales a ppear ^200-300 K cooler than 
would be inferred from theoretical model atmospheres (e.g. I Jones et al.l 120051 ). for example. 



While M dwarf masses, temperatures, and radii have proven challenging to calibrate, 
those parameters appear strikin gly pedestrian when contrasted with the difficulty of 



inferring M dwarf metallicities. 



Mould! ( Il976a[ ) was the first to systematically assess the 



strengths of M dwarf spectral features as a function of metallicity, constructing a grid of 
synthetic model atmospheres that included molecular opacities and spanned temperatures 
of 4250-3000K and metallicities of [Fe/H]= -2.0 - 0.0 dex. With these models, Mould 
established several foundations of subsequent investigations of M dwarf metallicities, such 
as the metallicity sensitivity of TiO absorption, the gravity sensitivity of CaH absorptio n, 



and the rough metallicities of Galactic M subdwarfs ( [Fe/H]~ -1.0 dex; [Mould 



1976bJ). 



These studies enabled the subsequent photometric and spe ctroscopic identification of 



M subdwarfs (IStauffer fc Hartmann 



1986 



Ruiz fc Anguita 



increasingly detailed models of M dwarf atmospheres ( lAllard 



1993), 



;he development of 



1990), and represent the 
fundament al origin of t he TiO and CaH-based indices used today to identify metal-poor M 



subdwarfs (IGizis 



1997 



Woolf & Wallerstein 



2005 



Lepine et al. 



20071 ). 



While techniques to identify metal-poor M dwarfs have been developing for more 
than 35 years, methods for identifying metal-rich M dwarfs have received significantly 
less attention. Interest in identifying M dwarfs with super-solar metallicities was only 
recently spurred by the realization that such stars could be highly promising targets 
for searches for Earth-like planets. The existence of a correlation between (gas giant) 
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planet frequency and host-star metallicity is we 



1997 



Santos ct al. 



2004 



Valenti fc Fischer 



1 established for FGK-dwarfs ( I Gonzalez 



20051 ). but the first investigations to examine 



the planet-metalli c ity correlation ext ended to M dwarf hosts returned mixed results 



( iBonfils et al. 



2005a 



Bean et al. 



2006b|). Bonfils et al. (2005, B05 hereafter) analyzed V 
and K photometry for 20 wide M dwarf companions to FGK-dwarfs whose metallicities 
could be determined reliably via standard spectroscopic techniques. Assuming that both 
binary components inherit the same metallicity from their parent molecular cloud material, 
B05 assigned the metallicity measurements for each primary to the M dwarf secondary, 
and derived iso-metallicity contours in the My vs. V-K s color-magnitude plane. This 
photometric metallicity calibration suggested that nearby M dwarfs, including the planet 
hosts Gl 876 and Gl 436, were slightly metal-poor compared t o the mean metalli city of the 



Galactic disk. This conclusion was reinforced by the work of 



Bean et al. 



(120060 ). whose 



spectroscopic analysis indicated sub-solar metallicities for Gl 436, Gl 581, and Gl 876. 

The persistence of the planet-metallicity correlation into the M dwarf regime has been 
suggest ed, however, by the recent re-calibration of the photometric metallicity calibration 



by I Johnson fc Appd (J2009J). Using six M dwarfs with wide, metal-rich FGK companions to 
update the My vs. V-K metallicity contours, Johnson & Apps (2009, JA09 hereafter) found 
that the B05 calibration systematically underestimated the compositions of these metal-rich 
stars. JA09 used their revised calibration to estimate the metallicities of six M dwarfs with 
planetary mass companions, concluding that M dwarf planet hosts are indeed preferentially 
m etal rich, just like FGK h osts. This conclusion was supported by the subsequent work 



of ISchlaufman fc Laughlinl (j2010l ). who updated the B05 and JA09 calibrations by using 
theoretical models to inform the functional form of the M# vs. V-K relation, and by 
inferring the mean metallicity of M dwarfs in the Solar Neighborhood from a volume-limited 
and kinematically-matc hed sample of F & G dwarfs in the Geneva-Copenhagen Survey 
f lNordstrom et all 12004( 1 . 
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The photometric metallicity calibration developed by B05 and extended by JA09 and 
SL10 has proven to be a valuable resource, but its dependence on precise V magnitudes and 
trigonometric parallaxes limits its utility to M dwarfs in the immediate solar neighborhood, 
at least for the remainder of the pre-Gaia era. In a recent contribution, we presented an 
alternative , spectroscopic techniqu e for estimating the metallicities of near-solar metallicity 



M dwarfs (IRojas-Ayala et al. 



2010 



hereafter RA10). This technique is capable of providing 
metallicity estimates with an accuracy comparable to that of the photometric technique, 
and requires only a moderate resolution K band spectrum of the M dwarf target, providing 
a significantly lower observational burden for empirically estimating the metallicities of 
distant, near-solar or super-solar M dwarfs. An example of the utility of this technique is 
provided by a recent paper on M dwarf planet hosts identified in the most recent Kepler 
data release by Muirhead et al. (2011; submitted) 

While RA10 provided a concise introduction to the spectroscopic K band metallicity 
indicator, space limitations prevented a full exploration of the technique. In this paper, 
we provide a full analysis of the motivation, calibration, demonstration and application 
of the technique using spectra we obtained for 133 Solar Neighborhood M dwarfs. In 
Section 2, we describe the composition of our sample, and the acquisition and reduction 
of our spectroscopic observations. In Section 3, we describe our measurements of the new 
H20-K2 (a modified version of the H20-K index utilized by RA10), Na I and Ca I spectral 
features upon which the K band metallicity technique depends. In Section 4, we analyze 
model atmospheres to demonstrate the H20-K2 index's insensitivity to metallicity, and use 
the H20-K2 measurements from our sample to calibrate the index as a spectral type and 
T e ff indicator. In Section 5, we use the H20-K2 index as the basis for a revised K band 
metallicity calibration; by incorporating an additional M dwarf binary into our sample 
of metallicity calibrators, and adopting a modified functional form to better normalize 
the temperature dependence of the Na and Ca features, we derive an updated K band 
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[Fe/H]metallicity indicator, as well as a new calibration for overall metallicity ([M/H]), 
along with robust estimates of the uncertainty in each calibration. In Section 6 we perform 
several sanity checks of the K band metallicity technique, demonstrating that the K 
band metallicity estimates preserve expected correlations between metallicity and Galactic 
kinematics as well as signatures of chromospheric activity. We conclude in Section 7 by 
using the K band metallicity estimates, in combination with other stellar parameters such 
as rotational velocity (v sin i), magnetic activity (as diagnosed by H a emission) and optical 
brightness, to identify a sample of nearby M dwarfs with particular promise for exoplanet 
surveys. We summarize our findings in Section 8, and present our full K band spectral atlas 
in an appendix (we also make our spectra available to the community online). 



2. Sample Selection Sz Observations 

We observed 133 nearby M dwarfs with declinations higher than -30°. Our sample 
consists of: 



Eighteen M dwarfs with wide (> 5" separation), common-proper-motion solar-type 



companions to serve as metal 
metallicity measurements by 



icity calibrators. The FG K-primaries have spectroscopic 



Valenti fc Fischer! (12005 



SPOCS Catalogue), obtained 



by fitting synthetic atmospheric spectra to their high-resolution, high signal-to-noise 
echelle spectra. The SPOCS [Fe/H] and [M/H] values for the FGK-primaries have 
been assumed to also describe the metallicities of their M dwarf companions. This 
assumption is justified if both binary components formed together, from the same 
well-mixed molecular cloud, and no mass transfer or dredge-up has occurred in the 
system. This assumption is supported by the measurements of the metallicities 



of binaries with two FGK components: 



Desidera et al. 



( 12004 



20061 ) find typical 



metallicity differences of <0.02 dex between components of binary systems. The 
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binary sy stems were se 
stars, the Poveda et al- 



erted from the iGliese &: Jahreisd (Il99lh catalogue of nearby 



(I1994T ) catalogue of n earby wide binary and mu ltiple systems, 



and the list of new HIPPARCOS binaries bv iGould fc Chanamel (120041 ) . The SPOCS 
metallicity values are shown in Table |5j The stars in the calibration sample have 
[Fe/H] metallicities from -0.69 dex to +0.31 dex, [M/H] metallicities from -0.49 dex 
to +0.25 dex, and spectral types from Ml to M6. 

One hundred and fifteen M dwarfs within 33 pars ecs, including 11 M dwarf planet 



hosts observabl e from Palomar M ountain: Gl 876 (IMarcy et al. 



mm , Gl 436 
200ffl. Gl 176 



(Johnson et al. 



(lApps et al. 



Butler et al 



Forveille et al. 



2004h. Gl 581 (IBonfils et al 



2005b). G 



1998 



De 



20091 ). GJ 1214 ( Oharbonneau et al. 



fosse et al. 



2009), Gl 649 



2O10h . HIP 57050 (aka G 



2010h . Gl 179 (Howard et al. 



1148 , lHaghighipour et al 



20101), and HIP 12961 (IForveille et al.lboill ). 



849 (IButler et al. 



2010), HIP 79431 



More than half are members of t 



le Northern 8 parsec 



selected from the LPSM catalog (jLepine fc Shara 



2005) 



sample; the remainder were 



Near-infrared spectra of these stars were obtainec 



on the Palomar 200-inch Hale Telescope (jHerter et al. 



with the TripleSpec spectrograph 



20081 ) during several observing runs 



between 2007 and 2010. TripleSpec at Palomar has no moving parts and simultaneously 
acquires 5 cross-dispersed orders covering 1.0-2.4 //m at a resolution of A/AA~2700. 

The observations were carried out using two different methods: 

1) The star was placed at 2 different positions along the slit, A and B. Four exposures 
were taken with an ABBA slit-nodding pattern. Sky-subtraction was performed by 
differencing A and B exposures in each paired nod. 

2) The star was placed at 5 different positions along the slit. This type of observation 
pattern was preferred since the array showed a significant number of bad pixels. Each of 
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the 5 exposures were sky- subtracted using a sky- frame made by median combining the 5 
exposures of the star. 

These two types of sky-subtraction methods also remove the emission from the 
instrument and the dark signal from the detector. 

Calibration data were taken at the beginning of each night to correct for the 
non-uniform response between detector pixels. Several exposures of 30 seconds were taken 
with the dome dark (lamp-off frames), and while illuminating the instrument with a white 
light source (lamp-on frames). A master-flatfield was obtained by subtracting the median 
of the lamp-off frames from the median of the lamp-on frames. 

The spectra were reduced using a modified version of the TEDI reduction software 



(IMuirhead et al. 



20 1 ih l 1 !. Each sky-subtracted exposure was divi ded by the normalized 



master-flatfield, wavelength calibrated and optimally extracted ( IHorne 



19861 ). The IDL 



pipeline walks through the stellar spectrum in each order of the calibrated star image, finds 
the maximum every 10th pixel, and fits a Gaussian to the slit profile. The Gaussian fit 
therefore contains all of the stellar signal. The spectral orders of TripleSpec are substantially 
curved, so the pipeline interpolates onto a rectilinear grid with a thin-plate-spline to correct 
for the slit tilt, using the spatial and wavelength solutions for the detector by Terry HerteJ^. 

The 1-D M dwarf spectra were telluric c orrected using obs ervations of an AOV star 



with the IDL-based code xtellcor _general by 



Vacca et al. 



(120031 ). The spectrum of an AOV 



star is almost free of metal lines: its observed near-infrared spectrum consists mainly of a 
featureless and smooth continuum superposed by H I absorption from the star's atmosphere, 
and telluric absorption features contributed by the Earth's atmosphere. A telluric spectrum 



1 http : / / www. astro . cor nell .edu / ~muirhead / ^Downloads 
2 herter@astro. cornell.edu 
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is obtained by removing the hydrogen lines of the AOV star using a high-resolution model 
of Vega. The target spectrum is then divided by the telluric spectrum constructed from 
observations of a AOV star obtained near in time and close in airmass to the target. 

Finally, the telluric-corrected spectra were flux-calibrated using their 2MASS K band 
photometry. The flux density of the target star is proportional to the data number count 
(Da) after the telluric calibration: 



C D, 



Then, the star's spectrum can be flux calibrated from its m^ magnitude, the K s 
spectral response function Rk s , and the K s -band average flux density of a m^ s =0 star F° Ks 
by finding C : 



C = / R KsdA F 10 _ . 4mKs 



(2) 



The K s magn itudes for all the 



( Skruts kie et al. 



targ ets were obtained from the 2MASS All-Sky Catalog 



2006( 1. The 2MASS K s spectral response curve and F^- 



of Point Sources^ 
are available online 

Table [1] presents the properties of each star in the sample, including the star's distance, 
V and K s magnitudes, the date the spectrum was observed, and the average signal-to-noise 
ratio (SNR) obtained in the K band continuum. The full sample of K band spectra are 
presented in Figures [TH2"Tt the sample is ordered in Table [1] and in Figures [111271 in order of 



3 http:/ /vizier, u-strasbg.fr/viz-bin/ VizieR?-source=II/246 

4 http:/ /www. ipac.caltech.edu/2mass/releases/allsky/doc/sec6 _4a.html 



increasing H2O-K2 index (see Section [3T2]) . Four stars appear to have sub-optimal telluric 
corrections (Gl 644AB, 829AB, 809, & 908); these stars are flagged in Tableland their 
spectra are presented separately in Figure [2HJ 



3. Spectroscopic Analysis 

3.1. Na & Ca EWs 

RA10 demonstrated that the 2.205 /im Na I and 2.263 fim Ca I lines can be used to 
estimate the metallicities of M dwarf stars. We measured the equivalent widths (EWs) of 
each of these lines for every star in our sample. The standard definition of the EW of a line 
is given by the following equation 



EW> 



/•A2 


I F(X) 1 




L F C (X)\ 



dX 



(3) 



Here, F(A) represents the flux across the wavelength range of the line (A2-A1), and 
F C (A) represents the estimated continuum flux on either side of the absorption feature. 

The EWs of the Na I and Ca I features were calculated with an IDL pipeline using the 
following approximation: 



EW A ~ 

8=0 

Equation H] is the Riemann sum expression of the integral in Equation [31 where the 
F(Aj) and F c (Aj) are the line flux and the estimated continuum flux of the wavelength 
interval AAj, respectively, and n is the number of intervals. The integration limits adopted 
for the Na I doublet and Ca I triplet features are shown in Table [2j 



AA,; 



(4) 
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Since the continuum of the M dwarf spectra is affected greatly by molecular /broad 
absorption features, a pseudo-continuum is calculated instead using regions adjacent to 
the feature of interest and free of any other atomic features. The pseudo-continuum for 
each feature is estimated from a linear fit to the median flux within 0.003 /im (~ 7 pixels) 
wide regions centered on the continuum points listed in Table [2j the median flux at each 
continuum point is also weighted by its unce rtainty in perform i ng th e fit. The continuum 



Forster Schreiberl (120001 ) for use in modest 



points in Table |2] were originally defined by 
resolution spectra of late-type giants. 

To estimate the uncertainties on the EW measurements, a Monte-Carlo approach was 
used. An IDL pipeline added random Gaussian noise (based on the signal-to-noise of the 
spectrum) to the star's spectrum and then calculated the EW of each feature using Equation 
HI The entire procedure was repeated 1000 times and the adopted EW uncertainties of the 
Na I and Ca I features for each star correspond to the standard deviations of the 1000 EW 
measurements of each feature. 

The values of the Na I and Ca I EWs measured for each star, along with their errors, 
are presented in Table [3j 



3.2. H 2 0-K2 Index 



In RA10, the H 2 0-K index by 



Covey et al. 



( I2010l ) was used to account for the influence 



of temperat ure on the st r ength s of the Na and Ca lines being used to diagnose each star's 



metallicity. 



Covey et al. 



( 120101 ) adopted this index to characterize the spectral types of 
highly reddened young stars from moderate S/N spectra, and optimized the H2O-K index 
to sample the redder portion of the K band and ignored weak atomic features (Mg I, Ti I) 



that were not typically visible in their spectra. Stars in our sample, by contrast, possess 
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low extinctions and are relatively bright, such that two of the regions used to calculate the 
H 2 0-K index show obvious atomic absorption features in our high-S/N spectra of early 
M dwarfs. To ensure our water-index measurements are not affected or biased by these 
features, we developed a modified H^O-index using two new regions that do not show any 
noticeable atomic lines. This H2O-K2 water index is defined as: 



(^(2.070 - 2.090))/(^(2.235 - 2.255)) 
2 (^(2.235 - 2.255))/(J r (2.360- 2.380)) 1 ' 

where (J 7 (a — b)) denotes the median flux level in the wavelength range defined by a and 
b, in microns. This new index represents the change in the overall shape of the spectra of 
M dwarfs due to water absorption from 2.07 /im to 2.38 /zm, with smaller values of the 
H2O-K2 index correspondin g to greater amount s of H2O opacity. The difference between 



the H 2 0-K index defined by 



Covey et al. 



(120101 ) and the H2O-K2 index is shown in Figure 



Uncertainties in the H 2 0-K2 index measurements were computed using the same 
Monte-Carlo approach used to estimate uncertainties in the Na I and Ca I lines; each 
H2O-K2 error estimate represents the standard deviation of 1000 H2O-K2 measurements 
after adding synthetic noise to the spectrum consistent with the S/N value of the observed 
spectrum. 

The values of the H 2 0-K2 index measured for each star are included in Table [3] along 
with their uncertainties. 



3.3. Expectations from theoretical models 

The metallicity calibration presented in RA10 was empirically derived, with the 
temperature and metallicity dependences of the Na I , Ca I , and H 2 features inferred 
largely from visual inspection of the calibrator stars. While RA10 used PHOENIX model 



-14- 



atmospheres to qualitatively demonstrate the metallicity sensitivity of Na I and Ca I , 
and the reduced metallicity sensitivity of the H 2 feature, no detailed exploration of the 
behavior of the Na I , Ca I , and H 2 features as functions of temperature or metallicity 
were performed. 

The lack of a quantitative exploration of the line strengths predicted by model 
atmospheres is partially due to the fact that, until recently, synthetic models struggled 
to matc h the infrared spectrum of M dwarfs, which is predominately dominated by water 



opacity ( lAllard fc Hauschildt 



19951 ) . The discrepancy was believed to be due to incomplete 
water vapor lists; several water opacity profiles were used through the years, all of which 
over-predicted the K band opacity a nd resulted in a lac k of flux in the s ynthe tic model 



when compared to observed spectra (lAllard et al. 



20101). 



Asplund et al. 



(120091 ) recently 



derived a new estimate of the solar oxygen abundance, however, which is a factor of 2 
lower than previous estimates and seems to have solved t he wa ter opacity discrepancy. 



The new BT-Settl-2010 synthetic spectra by 



Allard et al 



(120101 ) incorporate the new solar 



abundances, along with updated molecular line lists, and provide a better match to the 
spectral distribution of M dwarfs across the near-infrared than previous models. The 
BT-Settl-2010 models are available onlinej; using this interface, we constructed a grid of 
BT-Settl-2010 synthetic spectra with effective temperatures from 2200K to 4100K, and 
overall metallicities [M/H] equal to +0.5, +0.3, +0.0, -0.5, and -1.0 dex. This effective 
temperature range covers the whole M dwarf sequence, and the model grid covers a range of 
metallicities larger than that spanned by our calibration stars, or the metallicities estimated 
with Equation [5] for any of our target stars. All the stars analyzed here are known to be 
dwarfs and none show any spectral features indicating otherwise, so a value of log g=5.0 
was selected to provide the best agreement with empirical M dwarf gravity determinations 



5 http: / / phoenix.ens-lyon.fr/simulator 
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( [Fernandez et al-lbooa IPemorv et al.ll2009l l. 



To provide a quantitative indication of the feature strengths predicted by the theoretical 
models, we degraded the resolution of each theoretical spectrum to match that of our 
observed sample and measured the EWs of the Na I doublet and Ca I triplet, as well as the 
H 2 0-K2 index, using the same routines applied to our observed spectra. 

Figure [30] shows the spectra of three stars in our sample, ranging in spectral type from 
Ml to M7, superimposed by BT-Settl-2010 synthetic spectra degraded to the resolution 
of the TripleSpec data. The overall shape of the observed spectra, dominated mostly by 
water absorption, are well matched by the BT-Settl-2010 models. The Na I and Ca I 
atomic features used by the metallicity calibration fits, however, are noticably weaker in the 
BT-Settl-2010 models than in the observed spectra: even the highest metallicity synthetic 
spectra cannot reproduce the strengths of the Na I doublet and Ca I triplet in the spectra of 
Gl 324 B and LHS 2090. The [M/H] =+0.3 and T eff = 4100K BT-Settl-2010 model matches 
the strengths of the CO bands and the neutral Ti and Al lines for HIP 12961 quite well, 
but it shows stronger Ca I features at ~1.97 /im, and no Mg I line at ~2.11 /im. These 
discrepancies could be explained by unc ertainties in the oscillato r strengths and opacities 



adopted for these lines, as suggested by 



Rajpurohit et al. 



( 120101 ). Another explanation 



could be the size of the wavelength interval (resolution) chosen to calculate the synthetic 
spectrum. If the lines in question are not well-sampled at the resolution used to calculate 
the synthetic spectrum, convolving the synthetic spectrum with a Gaussian to degrade the 
resolution to match TripleSpec could cause them to be artificially 'weakened'. 

Figure |3T] shows the H 2 0-K2 index and EWs of the Na I doublet and Ca I triplet 
measured from the BT-Settl-2010 models, after degrading their spectral resolution to match 
that of TripleSpec, as a function of effective temperature and color-coded by metallicity. 
The Ca I triplet behaves as expected from observed data: its strength increases as both 
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temperature and metallicity increase. However, the Ca I triplet disappears from the 
synthetic models at effective temperatures lower that ~2900K, even for the supersolar 
models. In our high S/N empirical spectra, the Ca I triplet remains visible even at spectral 
type M8. 

The water absorption features predicted by the models also agree reasonably well with 
expectations from empirical spectra, with water absorption increasing monotonically with 
decreasing temperature, in a manner independent of metallicity, down to T e g ~3000K. 
Below 3000K, however, the models show a sharp metallicity difference: for the subsolar 
metallicity models, the H 2 absorption remains almost constant from 3000K to 2200K, 
while water absorption in solar and supersolar models continues to increase with decreasing 
temperature. The stark difference in the H 2 opacity between solar and subsolar metallicity 
models suggests the discrepancy in H2O absorption may be more likely a computational 
artifact than an astrophysical difference. 

The most puzzling behavior is shown by the EWs for the Na I doublet. Observationally, 
the strength of the Na I doublet increases as temperature decrea ses, at least until effective 



temperatures equivalent to an M6V star (e.g. 



Cushing et al. 



20051 1 . Strong Na I doublets are 



also seen in the late-type stars studied in this work. The only synthetic spectra that show 
this trend, however, are the solar metallicity spectra. The most metal-poor ([M/H] = -1.0) 
synthetic spectra, by contrast, show the strongest Na I absorption near ~3000K, with the 
Na I EWs decreasing to both hotter and cooler temperatures; the moderately metal poor 
([M/H] = -0.5) models show a similar behavior but with a strong jump near T c g=2900K. 
The super-solar models, by contrast, possess the strongest Na I lines at warm temperatures, 
but the Na I lines in these models begin to weaken at T e fj < 3600K, ultimately predicting 
Na I line strengths that are anti- correlated with metallicity (in the solar /super-solar regime) 
below T eff =3000K. 
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The behavior of the Na I doublet in the solar and supersolar synthetic models is not 
reproduced in our empirical spectra. The metal-rich calibrator, Gl 376 B ([M/H] 5-poc.s 
• 0.1 I . [Fe/H] S p O (75=+0.20), is a late-type star (M6) and it exhibits one of the strongest 
Na I doublets of the 18 calibrators. Since low temperatures favor the formation of molecular 
species, the anomalous Na I behavior in the synthetic models could be an effect of the 
treatment of molecular formation at these temp eratures. The most recen t study of alkali 



Loddersl (I1999h 



Loddersl (119991 ) calculated 



chemistry in cool atmospheres was conducted by 
the mole fraction, the amount of an element divided by the total amount of the elements, 
of neutral Na, ionized Na, and the molecules NaCl, NaH, and NaOH as function of 



temperature at 1 bar total pressure 
computations performed by 



o r coo l dwarf atmospheres (3300K < T < 1500K). The 
Loddersl ( 119991 ) indicated that the mole fraction of ionized Na 
decreases towards cooler temperatures, while the mole fractions of NaCl and NaOH increase; 
the mole fracti ons of ne u tral N a and NaH remain almost constant over this temperature 



range. Indeed, iLodderj (ll999f ) found that neutral Na is the most abundant Na-based 
species over this temperature range, with the abundance peaking near ~2000K and then 
beginning to d ecline as Na2S condensation weakens the Na features in brown-dwarf spectra 



rtLodders 



2OO2L 



It is possible that these complex chemical networks could play a role in 



the anomalous behavior of the Na I feature in the synthetic models: perhaps the most 
metal-rich models favor the formation of Na-based molecules toward cooler temperatures, 
for example, such that the neutral Na responsible for the Na I 2.205 feature depletes more 
rapidly in metal-rich models, or perhaps Na2S condensation occurs at higher temperatures 
in the super-solar models. These suggestions are highly speculative, however, and based on 
the number of odd behaviors exhibited for models cooler than 3000K, we suspect the root 
cause is as likely to be computational as astrophysical in nature. 
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4. Calibrating The H 2 0-K2 index as a spectral type and T e ff indicator 



Water is a heteroatomic molecule and the appearance of the absorpti on band reflects the 



press ure and temperature sensitivity of the formation of these molecules ( [Merrill fc Ridgway 



19791 ) . Numerous investigators have leveraged this fact to use near-infra red water absorptio n 



bands to diagnose M dwarf effective temperatures and spectral types. 



Kleinmann fc Hall 



( 119861 ) noted that the depression at 4778 cm 1 ( 2.096 jttm) due to water absorpt i on wa s 



Jones et al. 



(1994) 



largest for dwarf stars, and increases with decreasing spectral type, 
placed this calibration on an absolute T e ff scale, deriving effective temperatures for 
stars with spectral types ranging from M2V to M9V by comparing the obser yed H2O 



bandh ead at 1.34 fim to laboratory data for the water absorption coefficient. 



Ali et al 



(119951 ) similarly created a T e ff indicator using an index sampling the K band H2O feature 



in low resolution spectra at ~2. 15 and ~2.20 fim, deriving a linear relationship between 
effective temp e rature from low-resolution spectra of MK-dwarf standards. More recently, 



McLean et al 



(120031 ) calculated linear relationships between stellar spectral type and the 



depth o 
dwarfs. 



water absorption ba nds present in R ~2000 near infrared spectra of M, L and T 



McLean et al. 



(120031 ) found that their K band water index, H2O-D, was constant 
for spectral types M0-M4, at which point the index decreased (due to increased H 2 
absorption) towards cooler spectral types. 

To calibrate the H2O-K2 index as a well-defined T e ff indicator, we provide in Table H] 
the values of the H 2 0-K2 indices measured from the BT-Settl-2010 models and originally 
presented in Figure I3T1 As Figure [31] demonstrates, the BT-Settl-2010 models predict that 
the H2O-K2 index is a monotonic function of temperature, but independent of metallicity, 
for stars in the 3000 K < T e ff < 4000 K temperature range. As noted earlier, the lower 
bound of this metallicity-independent behavior is defined by the clear difference in the 
behavior of the H 2 0-K2 index as a function of T e ff for metal-poor and solar/metal-rich 
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stars, respectively. While this low-temperature behavior may be a computational artifact, 
the metallicity dependent differences at T e j^s > 4000 K are much more likely to represent a 
true astrophysical behaviors. That is, while the H 2 opacity appears to 'saturate' for even 
very low metallicities, such that the most metal-poor stars in the model grid do not appear 
to possess any less H2O opacity than the most metal-rich stars of the same temperature, 
this is not true at the highest temperatures, near the onset of the H 2 feature in the NIR 
spectral sequence. Rather, the H2O index values tabulated in Table H] appears to indicate 
that the H2O feature appears at higher temperatures for more metal-rich stars. Care must 
therefore be exercised in inferring the temperature of stars with very weak H 2 features; in 
addition to the greater difficulty of ensuring a high S/N measurement of a weak feature, 
the BT-Settl-2010 models suggest that even if high quality H 2 0-K2 measurements can be 
obtained for stars with weak H 2 absorption, there will be large astrophysical uncertainties 
introduced into the resultant T e ff estimate due to the strong metallicity dependence of 
those features at relatively high T e ffS. The gravity dependence of the H 2 0-K2 index, as a 
function of model effective temperature and fixed metallicity, is shown in Figure [321 The 
gravity values of log g= 4.5 and 5.0 are the standard values adopted for M dwarf models. 
Within the metallicity range explored with the BT-Settl-2010 models, the sensitivity of 
the H 2 0-K2 index to surface gravity is negligible for T e jj->3000 K, and throughout the 
whole effective temperature range studied for the solar and super solar metallicity models. 
The largest discrepancies in H 2 0-K2 index due to surface gravity are found at lower 
temperatures in the [M/H]=-1.0 models. Therefore, according to the BT-Settl-2010 models, 
the H 2 0-K2 index values are metallicity and gravity insensitive for 3000K <T e ff< 3800K. 

Effective temperatures for the 133 stars were inferred from their H 2 0-K2 index 
using the solar ([M/H]=0.0) BT-Settl-2010 models. For each star, we interpolated the 
measured H 2 0-K2 value onto the solar metallicity model [H 2 0-K2,T e //] grid to get a T e ff 
estimate and error from its H 2 0-K2 uncertainty. The values used for the solar metallicity 



-20 - 



[H20-K2,T e //] grid are shown in Table HJ The effective temperatures derived should not 
be considered extremely accurate or reliable, especially for the late type M dwarfs in the 
sample, since the BT-Settl-2010 models of T e jy<3000K show H 2 absorption discrepancies 
due to [M/H] that are more likely to be a computational artifact that an astrophysical 
difference. Due to the finite resolution of the model grid, where models have only been 
calculated on a grid with a spacing of AT e ^=100K, systematic errors of ~100K are 
expected. 

The dominan t spec tral sequence for M dwarf stars was established two decades ago by 



Kirkpatrick et al. 



(11991 



hereafter KHM), who used red-optical spectra covering 6300-9000 
Ato define a set of primary and secondary M dwarf spectral standards. KHM derived 
their classifications from a least-squares minimization technique that depends on individual 
spectral features as well as the overall spectral slope. In the KHM technique, the low 
resolution spectrum ( 18 A) is only normalized using almost featureless region near 7500 
Abefore performing the least-squares minimization, so the overall shape of the spectra is 
conserved. 



Sixty eight of the stars in our sample have K HM spectral 



types 



Research Consortium On Nearby Stars (RECONS, iHenry et al. 



reported by the 



20061 ) . enabling us to 



calibrate the H2O-K2 index as a proxy for KHM spectral type. We exclude from this 
calibration sample 14 objects with potentially ambiguous H2O-K2 indices: eleven are 
members of binaries or multiple systems in very tight orbits, such that the final K band 
spectrum mixes emission from all the stars in the system; we exclude the remaining three 
stars on the basis of their sub-standard telluric/flux calibration. For the fifty four single 
stars in our sample with high-quality spectra and KHM spectral types, we performed a 
linear regression between the measured H2O-K2 indices and KHM spectral types using a 
Bayesian approach that takes into account the errors associated with each of the variables 



f lKellvl 



20071 ). According to 



Kirkpatrick et al. 



fll991f ). the KHM spectral types have ~ 0.5 



subtype errors; the errors associated with each H 2 0-K2 index are listed in Table [3J This 
regression analysis indicated that a star's KHM spectral type, expressed as a numerical M 
subtype, can be derived from its H 2 0-K2 indices as: 



M subype = A + B (H 2 0-K2) (6) 

A = 24.699 ± 0.930 

B = - 23.788 ± 1.067 

RMSE(M subtype) = 0.624 

The correlation between KHM spectral type, H 2 0-K2 index, and the residuals of the fit 
for our 54 spectral type calibrators, are shown in Figure |33J The residuals indicate that the 
H 2 0-K2 index can predict each star's spectral type with a typical accuracy of 0.6 sub-types. 
As Figure [33] shows, however, some stars possess KHM spectral types that differ from the 
type predicted from their H 2 0-K2 index by >1 subtype. A careful inspection of each of the 
fifty four objects rules out poor quality spectra as the cause of the outliers, but reveals that 
for the more extreme cases, the overall shape of their K band spectra disagrees with stars 
of the same KHM spectral type. These discrepancies could be explained by the influence of 
metallicity on the strengths of double-metal molecular bands, a topic we return to in §5.41 

We have used this relationship to infer the KHM spectral type for each star in our 
sample from its H 2 0-K2 index; these types are listed in Table El and were used to order 
the sample from earliest to latest spectral types. To provide high quality template spectra 
for subsequent spectral type measurements, we also combined all the spectra in our sample 
in bins according to their decimal H 2 0-K2 types. For each subtype bin, each spectrum 
was normalized by its mean flux between 2.16 and 2.2 /im. The spectra then were all 



-22 - 



cross-correlated with the spectrum of one of the stars in the subtype bin to make their 
spectral features coincide. Finally, the wavelengths offsets were applied to each spectrum. 
At least three stars, and as many as thirty, were then median combined to make the master 
template for each subtype, with the exception of the M9 template, for which LHS 2924 is 
the only prototype in our sample. The templates from MO to M9 are shown in Figure EU 

5. Updated and Extended K band Metallicity Calibrations 

5.1. A Revised Empirical [Fe/H] Calibration 

RA10 presented a linear equation for estimating an (early-to-mid-) M star's [Fe/H] 
value based on the strengths of its Na I and Ca I lines, and its temperature as diagnosed via 
its H2O-K index. We present here a revised formalism to describe this relationship, where 
the revision has been motivated by three factors: 

• The adoption of the revised H2O-K index, as introduced above; 

• The addition of a new M dwarf star, Gl 166C, to the metallicity calibrator sample 
(Table [5]) following the publication of RA10. 

• The adoption of a new functional form to estimate the star's metallicity from its K 
band spectral features. 

The adoption of a new functional form for the K band metallicity calibration is perhaps 
the most meaningful change, and deserves further explanation. The linear expression in 
RA10 produced unrealistically high [Fe/H] values for some of the more metal-rich M dwarfs 
in this sample (e.g., LHS 3799; [Fe/H] ~ +0.64 according to the RA10 calibration). While 
it is necessary to include the water index in the relationship to remove the temperature 
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dependence of the Na I and Ca I features, the assumption that water absorption should be 
included as an independent variable is not necessarily correct. Extensive experimentation 
with various functional forms revealed that a better [Fe/H] calibration fit could be obtained 
if water absorption is used as a "weight" for the Na I and Ca I strengths rather than as a 
fully independent variable. 

We conducted multivariate linear regressions on the Na I and Ca I EWs, weighted 
by each star's H2O-K2 index, for the 18 metallicity calibrators to identify the best fit 
relationship to predict each star's [Fe/H] value. The resulting calibration is: 



lT? /XJ1 a T3 ^a Iew n Ca I E w 

[Fe/H] = A + B i^=K2 + ° H^O=K2 (7) 

A = - 1.039 ± 0.170 

B = 0.092 ± 0.023 

C = 0.119 ± 0.033 

RMSE([Fe/H]) = 0.141 



For completeness, we summarize in Appendix |A] the Residual Mean Square (RMSp), 
the Root Mean Squared Error (RMSE) and adjusted square for the multiple correlation 
coefficient (R 2 ap ) statistics for evaluating the uncertainties and predictive power associated 
1 a functional fit to a g i ven se t of calibration data. These quantities were introduced 



wit 

by ISchlaufman fc Laughlinl (|2010l ) to compare the quality of the photometric [Fe/H] 
calibrations. The [Fe/H] fit in RA10 has a RMSp ({Fe/H]) =0.02 and an R 2 ap ([Fe/B])=0.63. 
The addition of another binary system and the new independent variables maintained or 
improved each of these statistics with respect to the RA10 calibration: the RA10 [Fe/H] 
calibration and the revision presented here possess equivalent RMSp values, but the new 
[Fe/H] fit has a higher R 2 ap ( J R^([Fe/H] J?A11 )=0.67). 
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For the [Fe/H] model in Equation [7J the Root-Mean-Squared Error of Validation 
(RMSEy, Appendix [AJ is equal to 0.161 dex. The RMSEy value is equivalent to a ~ 70% 
confidence interval or 1-a, and is considered a sensible estimate of average prediction error. 
The 95% confidence interval is equal to ± 0.344 dex. 

The [Fe/H] residuals versus the dependent variables of the fit in Equation [7] and the 
SPOCS [Fe/H] values are shown in Figure [351 The residuals show mostly random scatter, 
and no evident structure can be seen as a function of the weighted line strengths (i.e. left 
and middle panels), and the computed [Fe/H] from Equation [7] (right panel). The "outlier" 
in the bottom left of the [Fe/H] residuals vs the SPOCS [Fe/H] values plot is the lowest 
metallicity star in the calibration sample, Gl 611 B, however, and demonstrates the pressing 
need for an expanded sample of metallicity calibrators to establish if any systematic errors 
are present, and if so, remove them. 

The EWs of the Ca I triplet and the Na I doublet features, weighted by the H2O-K2 
index, are plotted in Figure 1361 for all the M dwarfs analyzed here. Iso-metallicity [Fe/H] 
contours calculated from Equation [7J are shown as dashed and dotted lines. In this plane, 
there is a clear distinction between the M dwarfs in the calibration sample with metal-rich 
and metal-poor FGK-dwarf companions. The M dwarf planet hosts also all have [Fe/H]>0.0 
dex according to this relationship, with the exception of Gl 581 and Gl 649, whose [Fe/H] 
are equal to -0.10 dex and -0.04 dex, respectively. The [Fe/H] values predicted by Equation 
[7J for the full sample of stars analyzed here are listed in Table [3j the most iron-rich star in 
the sample is the planet host HIP 79431(+0.46 dex) and the least iron-rich is the flare star 
V* V1513 Cyg (-0.64 dex). 
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5.2. A New K band Overall Metallicity ([M/H]) Calibration 



In addition to iron abundance estimates, the SPOCS catalogue also provides overall 
metallicity [M/H] values for t he FGK-companions of st ars in the calibration sample. In the 



fitting procedure described in 



Valenti fc Fischerl (120051 ) . the [M/H] value is an individual 



model parameter, rather than a quantity constructed from individual element abundances. 
We performed a linear regression as described above, assuming the same functional form as 
adopted for the iron metallicity estimate, to calibrate the relationship between the water 
weighted Na I doublet and Ca I triplet and [M/H]: 



[M/H] = A + B + c S*]™L (8) 

1/1 H 2 0-K2 H 2 0-K2 V ; 

A = -0.731 ± 0.120 

B = 0.066 ± 0.016 

C = 0.083 ± 0.023 

RMSEQM/H]) = 0.100 

The K band [M/H] fit has a RMSp([M/B])=0.010 and an adjusted R square 
R 2 ap {[M/R})= 0.67. For the [M/H] model in Equation EJ the RMSEy is equal to 0.111 dex 
(~ 70% confidence interval or 1-cr), and the 95% confidence interval is equal to ± 0.237. 
None of the previous photometric or spectroscopic works on M dwarf metallicities have 
estimated [M/H] values, therefore, a comparison to other methods cannot be performed. 
The [M/H] residuals versus the dependent variables of the fit in Equation [8] and the SPOCS 
[M/H] values are shown in Figure [371 The residuals show just random scatter: , as with the 
[Fe/H] residuals, no structure is evident in the residuals as function of the normalized Na 
and Ca line strengths or the computed [M/H] values. 

Figure [381 shows the same Ca/H 2 0-K2 versus Na/H 2 0-K2 plane as Figure |36| but the 
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dashed and dotted lines represent iso-overall-metallicity contours calculated from Equation 
|SJ The M dwarf planet hosts all have [M/H]>0.0 dex, with the exception of Gl 581 and 
Gl 649 (the same stars which were subsolar in [Fe/H]) with [M/H] values equal to -0.06 
dex and -0.02 dex, respectively. The most metal-rich star in the sample HIP 79431 has 
a [M/H]=+0.33 dex and the least metal-rich V* V1513 Cyg has [M/H]=-0.45 dex. The 
predicted [M/H] values from Equation [S] for the stars in our sample are shown in Table El 

As the BT-Settl-2010 theoretical model atmospheres are calibrated as a function 
of their overall metallicity, we are now able to directly compare the feature strengths 
measured from our model grid against the iso-metallicity contours we inferred from the 
empirical calibration sample above. Figure [39] shows the EW of Ca I vs the EW of Na I , 
both weighted by H 2 0-K2 that we measured from the BT-Settl-2010 model spectra with 
temperatures equal to or higher than 3000K. This plane splits the synthetic models by 
metallicity, in a similar way as it also splits the calibration stars for the metallicity fits. The 
[M/H] =-1.0 dex models lie in the bottom- left corner, and the [M/H] =+0.5 dex models lie 
upwards and to the right of the [M/H] =-1.0 dex models, consistent with the direction of the 
offsets predicted by the isometallicity contours. However, since the synthetic models predict 
weaker Na I doublet and Ca I triplet features than are observed in the empirical spectra, 
the models do not directly align with the empirical iso-metallicity contours (depicted by 
dashed lines) obtained with the K band [M/H] fit of Equation E The [M/H] =-1.0 dex 
models are located above the [M/H] =-0.6 dex empirical contour, while the [M/H]=+0.5 
dex model points are far below the [M/H]=0.0 dex empirical contour. The odd behavior 
of the Na I doublet features of the supersolar models mentioned above becomes evident as 
temperature decreases in Figure [39] Subsolar and solar models of temperatures between 
3700K and 3000K appear as parallel lines in this plane, but the supersolar models intersect 
the solar model at 3000K. 
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5.3. The K band [Fe/H] Method Compared with Other [Fe/H] Techniques 

5.3.1. Statistics of the K band [Fe/H] and Photometric Methods. 

We now seek to compare the quality of our [Fe/H] fit, and the resultant values it 
predicts for our calibration stars, to similar evaluati ons of the [Fe/H] predictions produced 



by the photometric calibrations. 



Neves et al. 



(12011 



hereafter Nil) performed a similar 



exercise, comparing the RMSp, R 2 apl and the dispersion around the mean (rms) values of 
the photometric metallicity scales usi ng a sample of 23 M dwar f s with solar-type companions 



of known metallicity, and refined the iSchlaufman fc Laughlinl (120101 ) calibration based on 
their results. 

Unfortunately, for the same reason that Nil did not include the RA10 metallicity 
predictions in their analysis, because our calibrator samples do not overlap well, it is difficult 
to directly compare the statistics of our fit to those computed by Nil. Nil analyzed 23 M 
dwarfs with reliable metallicity estimates, V magnitudes, and trigonometric parallaxes; of 
these, only 7 are members of our sample of M dwarfs with moderate resolution NIR spectra, 
providing minimal leverage for evaluating the accuracy of the various techniques with a 
common set of calibrators. This demonstrates an important point: the statistics calculated 
above, by definition, describe how well a given formalism can reproduce measurements 
of a particular calibration set. The meaning of those statistics, therefore, is inextricably 
linked to the particular calibration set from which they were computed: comparing those 
statistics amongst various metallicity calibrations requires that the statistics be computed 
with respect to a specific, and shared, set of calibrators. 

To make matters worse, the veracity of the calibration set is as important, if not more 
so, as the number of calibrators. For example, while the B05 functional fit has overall 
better statistics (as we show below), some of the stars in its calibration sample have poor 
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photometry values. The JA09, SL10, and Nil calibrations improved the B05 calibration 
by including only stars with high precision photometry /parallax measurements, excluding 
several of the stars in B05, and including stellar population and kinematic considerations. 
Moreover, we find differences up to 0.3 dex in the [Fe/H] values assigned for the M dwarfs in 
our calibration sample and the values used to derived the photometric scales. For example, 
B05, SL10 and Nil assign a [Fe/H] value of -0.15 dex to Gl 250 B, while the SPOCS [Fe/H] 
value we adopt for this star in our own calibration sample is +0.14 dexu- As the metallicity 
calibrations were derived with reference to these [Fe/H] values, the choice of which set 
of [Fe/H] values to adopt as a benchmark will necessarily bias the comparison towards 
the calibrations derived from those [Fe/H] values. As a result, a truly fair comparison of 
the photometric and spectroscopic techniques will only be possible once the calibration 
sample has been significantly enlarged, and the veracity of the [Fe/H] values adopted for 
the calibration sample has been well established. 

Given the limitations outlined above, we conclude that the fairest comparison of the 
spectroscopic and photometric techniques that can currently be made is to qualitatively 
compare the internal quality of each fit, as benchmarked by its native calibration sample, 
as well as to compare the size and character of the calibration sample from which each fit 
was derived. Table [6] lists the rms, RMSp and R 2 ap values for each calibration technique, 
along with the number of predictors p and sample size N for each model. The uncertainties 
were estimated using the bootstrapping method described in Nil. 

The number of calibrators and the dispersion around the mean value evaluate the 
quality of the calibration sample that each technique is based on. Large samples (large N) 



6 We note that as this paper went to press, lOnehag et al.l ( 120111 ) reported a third metal- 
licity for Gl 250B of -0.05 dex, based on an analysis of high-resolution (R ~50,000) high 
signal-to-noise J band spectrum of Gl 250 B itself. 
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that cover a large range of values (large dispersion) are preferred. According to Table El 
the B50 calibration sample is the largest (iV=48), and also the most diverse (rms=0.483), 
followed by our calibration (RAll), SL10 and Nil, respectively. The JA09 calibration 
sample has a very low rms value since its composed only by high [Fe/H] stars. 

The RMSp and R 2 ap statistics evaluate the internal quality of each functional fit; they 
describe how well the values predicted by each fit match the measured for the calibrator 
stars that were used to construct the fit. Small RMSp values (small residuals) and R 2 
values closer to 1 are preferred (see Appendix |Aj) . The JA09 fit has the smallest RMSp, 
but its R 2 ap has a large uncertainty. Nil noted that the R 2 p is a noisy statistical estimator 
for small samples, however, and the JA09 fit only has 6 stars in its calibration sample. The 
BO 5 fit has the largest R 2 . 

In summary, the B05 calibration has the better statistics when its native calibration 
sample is used. However, the B05 calibration underestimates M dwarf metallicities, 
probably a result of using old visual magnitudes, instead of Johnson V-band magnitudes, 
and the lack of high [Fe/H] M dwarfs in their calibration sample. The JA09 fit has a very 
small calibration sample composed only of high metallicity M dwarfs, and therefore, not a 
reliable fit for [Fe/H] prediction for M dwarfs. We find that the RAll, SL10 and Nil fits 
have calibration samples of similar size, similar rms and RMSp values; our calibration has 
the highest R 2 p values of the three. However, as noted earlier, none of these statistics can 
evaluate the veracity of the input values of the calibration sample, which determines the 
accuracy of the predicted values by the functional fit, as the recalibration of the B05 fit 
demonstrates. Constructed an enlarged calibration sample, and demonstrating the veracity 
of the adopted [Fe/H] estimates, are key steps to enable a full and detailed comparison of 
the reliability of these techniques. 
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5.3.2. K band [Fe/H] Estimates Compared with Other [Fe/H] Techniques 



Supersolar metallicities for most of the M dwarf planet hosts are predicted 



K b and [Fe/H] fit as well as the p hotometric calibrations by 



and 



Johnson fc Appsl (120091 ) 



by t 



he 



Schlaufman fc Laughlinl ( 12010I ). While these three methods agree on the systematic 



metal- richness of planet hosts, the [Fe/H] values they predict do differ for each of the planet 
hosts (Table HO]) . 

Photometric [Fe/H] methods need accurate V magnitudes and distances, to calculate 
absolute K5 magnitudes, and thereby infer the metallicities of M dwarfs from the My 
vs V-K5 color- magnitude diagram. Forty eight of the M dwarfs in our sample have 



accurate Johnson-V magnitudes and parallaxes 



and Ks magnitudes by 2MASS ( jSkrutskie et al. 



by HIPPARCOS (Ivan Leeuwen 



20071 ). 



20061 ) . By including the M dwarfs with 



ohnson V magnitudes and parallaxes in the Yale Trigonometric Parallaxes Catalogue 



( Ivan Altena et al. 



20011 ). we can assemble a sample of ninety stars with the necessary data 



to compare the results from the new K band [Fe/H] fit and photometric [Fe/H] relations. 

Table [7] lists the ninety M dwarfs with Johnson-V magnitudes and parallaxes, with 
their [Fe/H] values obtained from the ph o tome tric calibrations and the K band [Fe/H] 



fit in Equation [7J Since the 



Neves et al. 



(12011 



refinement of the ISchlaufman fc Laughlinl (120101 ) calibration, we only included the later in 



photometric calibration is a marginal 



this comparison. Eight-six of the stars in Table [7J are also incl uded in the Palo mar-Michigan 



State University Nearby Star Spectroscopic Survey (PMSU; 



Reid et al 



19971 ). and have 



measu rements of the TiO and CaH feature strengths from their optical spectra. 



Woolf et al. 



( 120091 ) created a spectroscopic optical technique based on these indic es to estimat e [Fe/ H] 



Woolf et al. 



(12009!) 



metallicities for early M dwarfs. The [Fe/H] values predicted by the 
technique from the PMSU indices are also listed in Table [7J The value of [Fe/H] is listed 
for a star in Table [7] only if the parameters to calculate it satisfy the conditions of each 
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technique: 



Bonfils et al.l ( I2005al . hereafter B05): The [Fe/H] photometric relation by B05 is valid 
for stars with 4 < M K < 7.5, 2.5 < V-K s < 6, and -1.5 < [Fe/H] < +0.2. 



Johnson fc Appd ( 120091 . hereafter JA09): The [Fe/H] photometric relation by JA09 is 
valid for stars with 3.9 < V-K s < 6.6 and [Fe/H] > -0.05 dex. However, we decreased 
the [Fe/H] lower limit to -0.12 dex, which corresponds to the extrapolation done for 
the star Gl 832 in JA09. 



Schlaufman &: Laughlinl ( 12010I . hereafter SL10): The [Fe/H] photometric relation by 
SL10 does not explicitly mention any conditions. However, we only estimated [Fe/H] 
values for that stars with 3 < V-K s < 7, the color range in the color-magnitude 
diagrams in SL10. 



Woolf et al 



(120091 . hereafter W09): The [Fe/H] relation based on Ti O and CaH 



Lepine et al. 



molec ular absorption by W09 depends on the (no/Can parameter by 
( 120071 ). The W09 calibration is valid for stars with -1.5 < [Fe/H]< +0.05 and 
3500K < T c fj < 4000K. We used the stars with spectral types earlier than M3. The 
(no/CaH parameter categorizes the main sequence M stars into 4 classes ordered by 
decreasing metallicity: dwarfs (dM), subdwarfs (sdM), extreme subdwarfs (esdM), 
and ultrasubdwarfs (usdM). The metallicity class for each star is also listed in Table 
7. 



Figure 1401 shows the difference between the [Fe/H] values obtained by Equation [TJ and 
the metallicity methods listed before. 

The B05 and W09 relations assign lower [Fe/H] values to the stars with solar and 
supersolar metallicities according to the K band [Fe/H] calibration. The stars with [Fe/H] 
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differences within 0.15 dex have K band [Fe/H] metallicities between ~ +0.1 dex and ~ -0.5 
dex, and correspond to 50% and 58% of the stars in the B05 and W07 samples, respectively. 

The JA09 relation assigns higher [Fe/H] values to the metal-poor stars in the K band 
[Fe/H] calibration. The stars with 0.15 dex [Fe/H] differences (i.e. ~la of the K band 
calibration) have K band [Fe/H] metallicities between -0.3 dex to +0.4 dex, and represent 
the 55% of the stars with JA09 metallicities. 

The SL10 relation [Fe/H] values agree within 0.15 dex with the K band [Fe/H] 
calibration for a 56% of the stars with SL10 metallicities. However, the dispersion for 
metal-poor stars is larger than for the stars with high K band [Fe/H] values, when compared 
to the SL10 technique. The SL10 technique "underestimates" the values of some of the 
solar and supersolar K band [Fe/H] stars, but for low K band [Fe/H] values, it can either 
underestimate or overestimate that value by up to 0.6 dex. GJ 1116 A (M6) has a K band 
metallicity of [Fe/Hj^n = -0.12 dex, but the SL10 method predicts a value of [Fe/H]s£i 
= +0.89 dex. The SL10 [Fe/H] value for GJ 1116 A is very unlikely and it may be a result 
of bad V photometry or a bad parallax measurement. 

All the plots in Figure HO] show a similar slope, except the [Fe/H] SL10 - [Fe/Hj^n 
vs. [Fe/H]^4ii plot. The W09 relation was calibrated using the CriO/CaH parameter and 
metallicities inferred from Fe I an d Ti I abundances of K- dw arfs and early M dwarfs 



from their high resolution spectra ( I Woolf fe Wallerstein 



20051 ) . B05 acknowledge in their 



conclusion that their [Fe/H] values are consistent (whe re the paramenter spaces overlap) 



with the [Fe/H] values from 



Woolf fe Wallerstein 



( 120051 ). Since B05 and W09 results agree, 



it is not surprising that both plots look quite similar. The JA09 calibration depends on the 
same observables as the B05 calibration, though supersolar metallicity stars were used as 
[Fe/H] calibrators for the fit. This could explain why the JA09 plot has the same slope as 
the B05 and W09, but with an offset towards higher metallicities. The observables in B05 
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were also used in the SL10 calibration, however, the kinematics of the FGK-calibrators and 
theoretical models were considered in the SL10 calibration. For the supersolar metallicity 
stars, the SL10 shows just an offset towards lower metallicities when compared to the RAll 
calibration. For the metal-poor stars, the scatter is larger but there is not a distinct slope 
between SLfO and RAf f . 



5.4. Inferring potential optical metallicity indicators 



While M dwarfs are bright in the NIR, historically the optical regime has received more 
scrutiny as a means of studying M dwarfs. The molecules CaH, TiO, CaOH, and VO cover 
the spectral ran ge between 6300-9000 A. 



dependent (e.g. 



Woolf fc Wallerstein 



aO absorption bands are known to be metallicity 
20061 ). Titanium oxide bands grow stronger with 
increasing metallicity as well as decreasing temperature. Hydride absorption bands also 
increase in strength with decreasing temperature, but, since the number density of neutral 
hydrogen increases as metallicity decreases, they grow stronger with decreasing abundance, 
too. Therefore, hydride absorption bands dominate the spectra of metal-poor stars, and the 
ratio between TiO and CaH has be en used as a discriminant between disk and halo M stars 
(e.g. the subdwarf classification by iGizisI 119971 ) . 



One hundred and nine M dwarfs studied in our sample ar e part of the Palomar /Michigan 



State University Nearby Star Spectroscopic Survey (PMSU, 



Reid et al 



1995 



Hawley et al 



19961 ). The PMSU survey obtaine d optical spectroscopy of most of the M stars in The 



Third Catalogue of Nearby Stars (jGliese fc Jahreiss 



19911 1 . allowing spectral types to be 



estimated from TiO band-strengths, chromospherically active stars to be identified by H a 
emission, and metal-poor stars to be identified from the CaH strength relative to TiO 
strength. The overlap between our sample and the PMSU catalog allow us to test if a M 
dwarf's metallicity influences the spectral type obtained by the KHM technique. 
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Reid et al.l (119951 ) defined nine narrow band spectroscopic indices (5 TiO indices, 3 CaH 
indices, and 1 CaOH index) to measure molecular features in the 6200-7200 A wavelength 



range. Assuming that TiO band-strengths are primarily temperature dependent, 



ReidetaL 



( 119951 ) measured the full strength of the TiO band at ~7050 A (the Ti05 index) for 
a large nu mber of sta rs with K HM spectra l type s to calibrate TiO as a spectral type 



indicator. 



Gizisl (119971 1 used the 



Reid et al. 



(119951 ) index definitions to calibrate CaH2 and 



CaH3 indices as spectral type indicators, too. We calculate c 



Reid et al 



spec tral t ypes 



19951 ) and 



or the stars 



Gizisl ( 119971 ) fit 



in our sample that are also in PMSU using the 
relations. The spectral types derived from the Ti05 and CaH features agree for most of 
the sample. The difference in spectral type predicted by the optical indices and H 2 0-K2 
index as function of overall metallicity is shown in Figure |4TJ The Ti05 index predicts later 
subtypes than the H 2 0-K2 index for the solar and metal-rich stars, and earlier subtypes for 
the metal-poor stars. The CaH2 index predicts later subtypes than the H 2 0-K2 index for 
most of the stars. The CaH3 index tends to give earlier subtypes than the H 2 0-K2 index. 
However, considering their uncertainties, the spectral types derived from the TiO and CaH 
indices agree with the H 2 0-K2 index based spectral types. 

The values of Ti02, Ti03, and CaHl indices versus the H 2 0-K2 index of the stars 
in the PMSU survey are shown in Figure H2J There is a strong correlation between the 
water index and the TiO indices up to H 2 0-K2 ~ 0.8. For spectral types later than M6, 
VO decreases the strength of the TiO features at ~7050 A, and therefore, reverses the 
relation. The CaHl index saturates at Ti05 ~0.49 and , henc e, wa s not used a s a spectral 



type i ndicator for the optical calibrations by 



Reid et al 



jl995j ) and 



Gizisl ( 119971 ). 



Reid et al. 



( 119951 ) noticed that even when CaHl was weaker that the CaH3 index, it appeared to be 
more sensitive to "(presumably) abundance variations amongst the earlier type M dwarfs" . 
Figure H2 has M dwarfs with [M/H] > 0.0 dex and [M/H]< 0.0 dex color-coded with red 
and blue, respectively. The considerable dispersion in the CaHl vs H 2 0-K2 relation can 
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be explained by the metallicities of the stars. The CaHl index saturates at ~M3, due to 
the increasing TiO absorption in that region of the optical spectrum. Since both molecules, 
CaH and TiO, increase in strength towards lower temperatures, a ratio between these two 
molecules could be what CaHl is really measuring. If a given star has low metallicity, the 
TiO absorption is not going to be as significant as the CaH, for any given temperature. The 
CaHl feature will remain strong for low metallicities, which will correspond to low CaHl 
index values. But, if the star is metal-rich, the pseudo-continuum used to estimate the 
CaHl index is going to be affected by the strong TiO absorption. The CaHl feature will 
become less prominent, corresponding to high CaHl index values. There is a similar split 
between the metal-rich and metal-poor M dwarfs for the Ti02 and Ti03 relations, but it 
is less significant. The CaHl feature could be used as a metallicity discriminator for early 
solar-metallicity M dwarfs in the optical. 



6. Independent Tests of the K band [Fe/H] Fit: M-M Binaries, UVW Space 

Motions and Ha Activity 

6.1. Do M-M Binaries Return Similar Metallicities? 

The K band [Fe/H] calibration and the photometric methods rely on the assumption 
that stars in binary systems share the same metallicity, as a result from being formed 
from the same interstellar material and at the same time. This assumption, which allows 
an M dwarf secondary's metallicity to be inferred from measurement s of its FGK primary 



2004 . 



has b een validated for binaries with multiple solar type components (iDesidera et al.l 
20061 ). and implies that any metallicity estimation technique should return consistent values 
for two M dwarfs in a wide multiple system. As a sanity check, therefore, we inspect the 
[Fe/H] values predicted by each technique for the individual components of double M dwarf 
wide binary systems in our sample. Table [8] shows the K band and photometric metallicities 
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of 5 M-M binary or multiple systems in our sample, along with their NIR spectral types; 
Gl 643 and Gl 644C are also included as a pair, since both are pr oper motion companions 



to the Gl 644AB pair/triplet according to 



Segransan et al 



( 120001 ) . Photometric estimates 



cannot be made for all the stars in the systems since their photometry and/or parallaxes 
do not satisfy the conditions listed in Section 15.31 for each calibration. We also show these 
systems in Figure HH1 where we connect the individual components of a multiple system 
with dotted lines to examine how well each system agrees with the slope of our derived 
iso-metallicity contours. 

The K band [Fe/H] and [M/H] estimates for individual components of (non-calibrator) 
binary systems agree remarkably well; four of the six systems have [Fe/H] and [M/H] 
metallicity differences of less than 0.02 dex, and none of the six have individual metallicity 
estimates that disagree by more than 0.13 dex, comparable to the la accuracy of the 
underlying calibration itself. Indeed, the largest discrepancy for the K band calibration 
is for the Gl 725AB pair, whose 0.09 dex and 0.13 dex [M/H] and [Fe/H] differences are 
similar to the dispersion (RMSE) of each fit in Equations M and [7J respectively. 

The metallicity estimates predicted by the photometric techniques, however, do 
not reproduce this agreement: the photometric techniques predict metallicities for these 
individual binary components that typically differ by 0.1-0.2 dex, comparable to the single 
largest metallicity difference predicted by the K band technique. The Gl 412 system is 
the binary in our sample with the greatest spectral type difference, consisting of Ml and 
M6 components. Despite the significant range in temperatures between these two stars, 
the metallicities predicted by the K band calibration only differ by 0.01 dex, suggesting 
our calibration is not strongly biased with respect to T e ff- The GJ1245 system provides 
another a good case study: both components have the same spectral-type in the optical 
and in K band, and their K band spectra share a very similar morphology, so it is perhaps 



-37- 



not surprising that the metallicity estimates predicted by the K band calibration for 
the individual components agree quite well (A RAll [Fe/H] = 0.01). The photometric 
[Fe/H] estimates for each component, however, disagree much more (A [Fe/H] = 0.25-0.3), 
reflecting underlying differences in the source photometry: while the distance measurements 
for GJ 1245 AC and GJ 1245 B agree well (5 d = 0.2 pc), their photometry does not, 
with GJ 1245 AC appearing ~0.55 mag brighter in both V and K s , placing it significantly 
higher in the My vs. V-K s diagram and thu s imply ing a more metal-rich composition via 



the photometric technique. 



McCarthy et al. 



(119881 ) resolved each component of the GJ 



1245 triple system with speckle inte r 



of its components. 



McCarthy et al. 



erom etry and obtained K band photometry for each 



( 119881 ) 's photometric measurements revealed almost 
identical magnitudes for the A and B components, ~1.1 mag brighter than C. However, the 
combined K magnitude for the GJ 1245 AC system was ~ 0.5 mag brighter than GJ 1245 
B, but the optical spectra for AC and B were almost identical. Therefore, the faint GJ 1245 
C component could explain the photometric [Fe/H] discrepancy: it makes the star GJ 1245 
A to appear brighter without affecting its optical and K band spectra. This implies that 
an unknown unresolved binary could appear to be more metal-rich under the photometric 



calibration if it is assumed to be a single star. 



6.2. Do K band Metallicities Reproduce the Metallicity /Kinematic 
Correlations Expected for Galactic Sub-populations? 



It has been known for the better part of a century that there are clear correlations 
between the spatial, kinematic, and chemical properties of stars in the Milky Way. These 
correlations allow the Milky Way to be de composed into d istinct stellar populations and 



structural components (e.g., halo vs. disk; 



O'Connell 



19581 ) . although some debate remains 



concerning the exact number and properties of sub- components (e.g., thin vs. thick disk; 
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Ivezic et al.ll2008l ). iBensby &: Feltzingl (120 lOf ) recently showed that when 899 F- and G- 
dwarf stars were separated by their kinematic properties, a considerable number of thick 
disk stars exhibit metallicities consistent with the metallicity trend of the thin disk, and vice 
versa. When the solar-type stars were separated by age rather than kinematics, however, 
the number of outliers in each metallicity trend was reduced. 

Age estimates for M dwarf stars are difficult to obtain, however, making kinematic 



properties the best available means 



components. iBochanski et al 



or associating low-mass stars with various Galactic 



( 120071 ) investigated metallicity differences between M 



dwarfs associated with the thi n and thick d i sk pop u 



using the method described by 



Bensby et al. 



(I2003h . 



the (CaH2+CaH3)/Ti05 ratio, previously shown by 



ations, bisected kinem atically 



Bochanski et al. 



Lepine et al 



^2007J) measured 



to roughly 



discriminate between solar-metallicity ([M/H]~0), subdwarf ([M/H]~-1.2) and extreme- 



sub dwar: 
spectra. 



f[M/H1~-2.0) M dwar fs, for 6577 M dwarfs with Sloan Digital Sky Survey (SDSS) 



Bochanski et al. 



(120071 1 found that the observed (CaH2+CaH3)/Ti05 distributions 
for the kinematically divided thin and thick disk stars did not differ greatly, implying 
metallicity differences less than 1 dex, but suggested that the kinematically selected thick 
disk stars were consistent with an older population due to the lower fraction of stars with 
H Q emission, a tracer of chromospheric activity. 

One hundred fourteen M dwarfs studied in this work have UVW space motions 
from the Paloma r-Michigan State University Nearby Star Spectroscopic Survey (PMSU; 



Reid et al 



19971 ). enabling us to test if the K band metallicity calibration can replicate 



the metallicity trends identified for higher-mass members of the kinematic sub-populations 
of the Milky Way. One way to visualize the different kinematic populations of the Milky 
Way is through the so-called Toomre diagra nfl The Toomre diagram is a representation of 



Sandage fc Fouts 



(119871 ) were the first to show this diagram in their paper on the collapse 
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the stars' combined vertical and radial kinetic energies, (U 2 + W 2 ) 1//2 , as a function of the 
stars' rotational energy, V. Constant values of the total space velocity for the stars, V to t= 
(U 2 + W 2 +V 2 ) 1//2 , can be represented by circles centered at [0,0], as is shown in steps 
of 50 km/s in Figure Ml W e base our kinematic sub-divisions on the boundaries defined 



by iBensby &: Feltzing! ()2010l ). where thin disk stars possess Vi O £<50 km/s and thick disk 
stars possess ~70 km/s <V to £<200 km/s; we additionally classify stars with intermediate 
velocities (50 km/s <Vi O £<70 km/s) as thin-thick stars. We show in Figure H5] the 
distributions of K band overall metallicity estimates for members of each of these kinematics 
groups, with the mean and median [M/H] indicated for each distribution. Consistent with 
results seen for higher mass stars, we find the thin disk population is noticably enriched (A 
[M/H] ~ 0.2 dex) relative to the thick disk population, and a two-sided K-S test identifies 
that a <0.003% probability that the two samples are consistent with the same parent 
distribution. 



To further explore the differences between these populations, we show in Figure 
the Toomre diagram for the stars with kinematic measurements in the PMSU survey, 
color-coded by metallicity, with boundaries for the metallicity sub-groups chosen from the 
mean values of the metallicity distributions of Figure [451 Visualizing the sample in this 
way, it is apparent that the metal-poor stars span a much larger range of velocities than the 
metal-rich stars: more than half the stars with [M/H] < -0.14 have V to t > 70 km/s, while 
every star with [M/ H] > -0.15 lies with i n the V^^lOO km/s contour. These results reinfo rce 



the tr ends seen by 



Bensby fc Feltzing) (120101 ) for higher mass stars: 



Bensby fc Feltzing 



(120101 ) found that older (t > 8 Gyrs) stars are a kinematically hot sample that cover most 



oft 



re velocity space in the T oomre diagram, just like the thick-disk M dwarfs in Figure 



Bensby fc Feltzing! (120101 ) also found kinematically hot stars within the younger (t < 7 



of the Galaxy, but acknowledged the diagram to A. Toomre [1980] 
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Gyrs) population, but the majority have Vtot< 100 k/m, consistent with what we see in 
Figure HH for our (metallicity selected) thin-disk M dwarfs. 



6.3. Are K band Metallicities Consistent With an Assumed 
Age-Metallicity- Activity Relation? 

A final observable that can be used to probe the existence of an underlying age- 
metallicity relation is chromospheric activity. Common tracers of chromospheric activity, 
such as emission in Ca II H & K or H a , indicate that young stars demonstrate high 
levels of chromospheric activity, w ith the activity d ecaying as the stars age; this effect 



has been seen in solar type stars (jSkumanich 



(jStauffer fc Hartmann 



19721 ) and well into the M dwarf regime 



19861 ). Chromospheric activity timescales of M dwarfs expand from 



lGyr for early- type /higher-mass st ars, to ~ 10 Gyr for the later-type/lower-mass stars 



(IHawley et al. 



1999 



West et al 



2008|). 



As noted in the previous section, iBensby fc Feltzingl ( 120101 ) found that the link between 
age and metallicity appears to be stronger for solar-type stars than the link between 
kinematics and metallicity. While it is not possible to derive an exact age for an M dwarf 
based on its observed activity level, one can use the act ivity lifetim e s not ed above to 



Gizis et al 



)2) 



(120021 ) measured the 



statistically divide stars into older and younger samples, 
equivalent widths of H a on high-resolution spectra of 676 M dwarfs, and identified distinct 
groups according to their H„ properties. Our sample includes eighty two stars analyzed by 



Gizis et al. 



( 120021 ) . providing an opportunity to use activity as a crude proxy for age, and 



examine the evidence for an age-metallicity relationship within the M dwarfs in our sample. 

Figure H6] shows the H a emission strength versus spectral type for these stars, divided 
according to the metallicity boundaries inferred from Figure H5] and applied previously in 
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Figure HH Across all metallicity categories, most of the early-type (MO to M2) stars show 
H Q absorption, implying either a weak/moderate activity or no activity at all, while the 
later-type stars show H a emission, consistent with the longer activity lifetimes expected for 
the lowest-mass stars. The metal-rich category seems to have a slightly higher average level 
of H Q emission when compared with the other categories, consistent with the metal-rich 
stars being a statistically younger sample than the more metal-poor samples. If we adopt 
a 2 A threshold for significant H Q activity, the earliest star in each category with elevated 
activity is clearly correlated with metallicity: the earliest active metal-poor star is an M5, 
while the earliest active intermediate and metal-rich stars are M4 and M3, respectively. This 
correlation is precisely what would be expected if there were an underlying age-metallicity 
relation, such that the sample of metal-poor stars consists of statistically older stars, where 
activity is only seen for the lower-mass stars with relatively long activity lifetimes, while 
metal-rich stars are statistically young, with higher mass stars not yet having transitioned 
into the inactive stage. 



7. Particularly promising M dwarf planet hosts 



The nearly 600 confirmed extrasolar planets discovered to date reveal that planets are a 
natural and frequent by-product of star formation. One of the most important observational 
constraints is that planet freq uency rises steeply with host star abundanc e above solar 



metallicity in FGK-stars (e.g. 



Santos et al. 



2004 



Fischer fc Valenti 



20051 ) . As mentioned 



in the introduction, the first photometric and spectroscopic metallicity measurements 



for M dwarf planet host s indicated that they were slightly metal- poor ( Bonfil s et a 



2005a 



Bean et al 



Rojas-Ayala et al. 



2010 



2006a ). However, recent meta 



Schlaufman fc Laughlin 



licity calibrations (IJohnson fc Apps 



2009 ; 



20101 ) have convincingly shown that the 



small sample of nearby M dwarf planet hosts have solar and super-solar metallicities. 
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Table [3] flags the M dwarf planet hosts in our sample. According to the metallicity 
calibration derived in Section 5.2, the eleven M dwarf planet hosts have overall metallicities 
higher than -0.06 dex, with the Jovian hosts being more metal-rich than the Neptune 
or super-Earth hosts. The sample of M dwarf planet hosts is small, but their K band 
metallicities suggest the existence of a planet-metallicity correlation for M dwarfs. We 
compared the [Fe/H] distribution of the eleven planet hosts to the [Fe/H] distribution of 
35 M dwarfs in our sample that are also members of the California Planet Survey (CPS) 
sample of low- mass stars (John Johnson, private communication). These 35 M dwarf have 
been subjected to an intensive RV monitoring campaign by the CPS team, and their RV 
measurements have ruled out the presence of Jupiter-size planets within several AU. Figure 
1471 shows the metallicity distributions of the 35 CPS M dwarfs and of the 11 planet hosts 
in our sample. A two-sided K-S test identifies that there is a <1.6% probability that 
the K band metallicities we calculate for these samples of M dwarfs with and without 
planets share the same parent distribution. If the planet sample is divided by Jupiter and 
Neptune/super-Earth hosts, then a two-sided K-S test gives a 0.8% chance that the Jovian 
hosts and the 35 CPS stars are drawn from the same parent distribution. The probability 
that the Neptune and super-Earth hosts share the same [Fe/H] distribution as the 35 CPS 
stars without detected planets is much higher: 8.6%. These results are in agreement with 
what is observed for FGK stars, where Neptunes and super-Earths, unlike Ju piters, do not 



seem to form preferentia lly around highly metal-rich ([Fe/H]> 0.2 dex) stars ( ISousa et al 



2008 



Mayor et al. 



20111 ) . However, we acknowledge that the small size of these samples 



prevents any conclusive result. 

According to the [Fe/H] distribution of the planet hosts in our sample, M dwarfs with 
solar and super-solar metallicities should be the preferred targets to look for planets around 
cool stars. The most metal-rich stars ([M/H]>+0.2 dex) in this sample are the Jovian hosts 
HIP 79431 and Gl 849, and the M dwarfs Gl 285 (M5), G 203-47 (M4), Gl 169.1 A (M3), 
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LHS 3799 (M5), and Gl 205 (MO). Bower et al.l (120091 ) used radio astrometry to search for 



planets around a sample of active M stars within 10 pc and their a nalyzed samp 



two of the metal-rich stars listed above (Gl 285 and LHS 3799). 



Bower et al 



e mc 



20091 ) 



uded 



were able to rule o u t the existence of 3-6 Mj planets within 1 AU at the 99% confidence 



level. 



Lawler et al. 



(120091 ) searched for an infrared excess around Gl 205, and found none, 



excluding the presence of debris at distances beyond the snow-line. 

Interestingly, the metal-rich stars G 203-47 and Gl 169.1 A do not have planetary 
companions, but they do possess white-dwarf (WD) companions. Gl 169.1 B is a relatively 
young and featureless white dwarf wit h an age estimate o f ~ 2.8 Gyr based on a model 



with a core composed of pure carbon 



age estimate, however 



Delfosse et al. 



(jLeggett et al. 



19981 ). Gl 203-47 B does not have an 



(119991 ) predicted that its primary (Gl 203-47) would 



have a peculiar composition since the system was likely in a contact configuration during 
the post main-sequence evolution of the current day WD component. The sem i-major axis 



19991 ). Therefore, 



of the present orbit of the G 203-47 system is ~ 0.05 AU (iDelfosse et al.l 
G 203-47 could have accreted nucleosynthesis products dredged up to the surface of Gl 
203-47B's AGB progenitor. Figure 1481 shows the optical PMSU spectra and K band spectra 
of G 203-47 and Gl 169.1 A, along with the spectra of LHS 3409, a star with similar spectral 
type but subsolar metallicity. The strong Na I and Ca I features in K band, and the weak 
CaH absorption in their optical spectra, indicates that these stars have indeed supersolar 
metallicities. 

Most of exoplanet searches have been performed at visible wavelengths. M dwarfs are 
intrinsically faint in the optical, so the largest stellar component in the solar neighborhood 
haven't yet been searched thoroughly for planets, except for the brightest, closest early 
M dwarfs. However, plan et searches are now being performed with near- infrared rad ial 



velocity instruments (e.g. 



Bean et al 



2010 



Crockett et al. 



2011 



Muirhead et al. 



20111 ) and 
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with far red filters for transits (jlrwin et al.ll2009l ). For these searches, the visible faintness of 
M dwarfs is not a problem, but other intrinsic properties of M dwarfs should be considered 
to optimize detectability. 

Therefore, the most "desirable" M dwarf targets for planet searches should be: 

• Bright: The optical faintness of M dwarfs is not a problem for the near-infrared 
searches, but a large number of M dwarfs per night is desirable. Brighter M star are 
favorable for achieving a fixed signal-to-noise ration and for ease of follow-up for any 
planets that are detected. 

• Slow rotators: Stars with high rotation rates have radial velocity me asurements 
with 



Bouchy et al. 



ower precision since stellar rotation widens the spectral features. 
(j200ll ) found that the radial velocity uncertainties increase by a factor of ~4 when 
the rotation velocity increases from ~1 km/s to 10 km/s. Stars with v sin i> 10 km/s 
should be avoided to provide optimal RV sensitivity. 

Inactiv e: The activit y of M dwarfs can also cause false positives (e.g. starspot-crossing 



events; 



Carter et al. 



(I Wright 



2005 



201 



Bertaetal 



)) an d can cause radial velocity and photometric jitter 



20111 ). Depending on the starspot cover age, the planet-to- 



star radius ratio can change up to few percent (ICarter et al 



20111). Stars with H c 



emission should be avoided to avoid false positives and jitter. 

Metal-rich: The M dwarf planet host shown in Table [TTJ1 all have metallicities [M/H]> 
-0.05 dex. Therefore, it is probable that solar and supersolar metallicity M dwarfs are 
more likely to host planets. 



Table IH] lists the stars that satisfy the conditions mentioned a 
and metallicity in our M dwarf sample. The H Q values are from 



Dove for ac 



Gizis et al. 



ti yity, r otation, 
J2OO2J). The 
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v sin lvalues were taken from 



Jenkins et al 



(|2009| , Tables 1 and 3) and 



Browning et al 



(120101 ) . Some of the stars in Tabled have detected stellar or planetary companions already 
(five of the eleven M dwarf planet hosts have their H a and v sin ^listed in the papers 
mentioned above). It will be interesting to see how many planets are discovered around the 
targets in Table [9] in the near future, and what "flavor" they might be. 



8. Conclusions 

We present K band spectra for a sample of 133 M dwarfs, including 18 M dwarfs with 
reliable metallicity estimates (as inferred from an FGK type companion), 11 planet hosts, 
more than 2/3 of the Northern 8 pc sample, and additional M dwarfs from the Lepine Shara 
Catalog. From these spectra, we measured EWs of the Ca I and Na I lines, as well as an 
index quantifying the absorption due to H 2 opacity, features we have previously shown to 
be of use for predicting M dwarf temperatures and metallicities. 

1. From the subset of our stars which possess primary or secondary KHM spectral types, 
we calibrated the H 2 0-K2 index as a spectral type indicator. We also estimated 
effective temperatures for the stars in the sample by interpolating their H 2 0-K2 
indices onto a solar metallicity [H 2 0-K2, T e ff] model grid. 

2. We revise the functional form adopted to predict the [Fe/H] metallicities of M dwarfs 
based on measurements of their Na I , Ca I , and H 2 features. We perform a 
linear regression upon 18 M dwarfs with reliable metallicity estimates from FGK 
companions to calibrate this updated relationship. Statistical tests demonstrate that 
the [Fe/H] estimates produced by this relationship are accurate to RMSE =0.141 dex, 
and confirm that this relation accounts for more of the variance within our calibrant 
sample than any other existing technique for estimating M dwarf metallicities. 
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3. For the first time, we derive an expression for an M dwarf's overall metallicity [M/H]; 
quantitative comparisons to model atmospheres, which are benchmarked according 
to overall metallicity, provide a qualitative validation of our approach for estimating 
metallicities from Na I , Ca I , and H 2 features, but also demonstrate an overall 
offset between the atomic line strengths predicted by models as compared to actual 
observations. 

4. We examined previous optical molecular indices sensitive to stellar metallicity using 
out metallicity estimates. We identify that the CaHl feature as a potentially valuable 
optical metallicity discriminator for solar-metallicity early type M dwarfs. 

5. We perform several sanity checks of our metallicity estimates, confirming that our 
metallicity estimates reproduce expected correlations between metallicity and Galactic 
space motions and Ha emission line strengths, and return statistically identical 
metallicities for members of M-M multiple systems. 

6. We recovered the results from previous metallicity studies that nearby M dwarf planet 
host exhibit solar to supersolar metallicities, where stars with Jovian-mass planets 
are more metal-rich than stars with Neptune-likes or super-Earths. A list of the 
best targets in the sample for planetary searches selected by metallicity, activity and 
rotation-rates is given. 
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A. Statistical tests and reliable uncertainties for functional fits 

A well-fitting regression model results in predicted values close to the observed data 
values. The following statistics are commonly used to evaluate a model fit: 

Residual Mean Square (RMSp) 

The RMSp is defined as 



RMSp = ( A1 ) 

(A2) 
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where n is the number of data points (in this case the number of calibrators), p is 
the number of predictors in the model (including the constant term), y^ is the value of the 
response variable i and % is the prediction value given by the regression model for y^. The 
RMSp is the variance of the residuals. In general, an RMSp value closer to indicates a 
model fit that is more useful for prediction. 

Root Mean Square Error (RMSE) 

The RMSE is the square root of the variance of the residuals 



The RMSE indicates how accurately the regression model predicts the response. The 
RMSE is an absolute measure of it and has the same units as the response variable. Lower 
values of RMSE indicate better fit. 

Adjusted Square of the Multiple Correlation Coefficient {R 2 ap ) 

The R 2 p represents the proportion of variability in a data set that is accounted for by 
a regression model. The R 2 is defined as 



RMSE = VRMSp 



(A3) 



(A4) 



ap 



1 - 



(n — 
RMSp 
MST 




p) EILi (yi - y) 2 



(A5) 



1 - 



(A6) 



where y is the overall mean of the observations yi, and MST is the total mean of 
squares. The R 2 ap compares the unbiased variance of the residuals (RMSp) and the unbiased 
variances of the observations (MST), therefore, can be interpreted as the proportion of total 
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variance that is explained by the model. A value of R 2 ap =l indicates that the regression 
model explains all of the variance in the sample, while R\ =0 indicates that the regression 
models explains none of the variance. R 2 ap should always be used with models with more 
than one predictor variable. 

Root Mean Square Error of Validation (RMSEy) 

A better estimate of the uncertainty in the K band [Fe/H] fit can be obtained by 



Weisbere 



2005|) 



performing the Predicted Residual Sum of Squares procedure (PRESS; 
The PRESS statistic is equivalent to "leave-one-out" method, where, in this case, the 
regression procedure is repeated leaving one of the n calibrators out and performing the 
regression with the other n-1 calibrators, for each one of the n calibrators. The PRESS 
statistics is defined as 



PRESS = J^ef (A7) 

i=l 
n 

i=l 

where ei is the residual for the calibrator Cj, computed as the difference between the 
observed value of the predictand yt and the prediction yr-a)i from a regression model 
calibrated using the calibration sample but without the calibrator Cj. Then, it is possible to 
use as a measure of uncertainty a statistic based on this validation data 



RMSE V = (A8) 
V n v 

Equation (2.6) is known as the Root-Mean-Squared Error of validation (RMSEy), 
where ny is the number of calibrators. The RMSEy value is a sensible estimate of average 



prediction error according to iWeisbera (120051 ). The RMSEy value can be used then to 
obtain confidence intervals at a desired significance level around the predictions 



% ± C(k /2 , n _p) RMSE V 

where the multiplier of RMSEy corresponds to the desired probability point on the 
cumulative distribution function. The correct cumulative distribution function for this 
specific fit is the two-sided t-student distribution with n-p=15 degrees of freedom, where 
n=18 is the sample size for calibration and p=3 is the number of regressors in the model fit, 
including the constant term. The two-sided i-student distribution is preferred instead of a 
normal distribution due to the small number of calibrators (n= iiy =18). For approximate 
70%, 95%, and 99% confidence intervals, the values for C(t a /2,n-p=i5) are equal to 1.07, 
2.13, and 2.95, respectively. 
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Fig. Set 1. K band spectra of 133 nearby M-dwarfs 
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Fig. 1. — TripleSpec 200-inch K band spectra of 133 nearby M-dwarfs. 
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Fig. 29.— The H 2 0-K index by 



Covey et al 



( 120101 ) versus the new H 2 0-K2 index in Equation 



[5j The portions of the spectrum used by the H 2 0-K2 index are almost free of atomic features 
and better sample the change in the overall shape of the spectra of M dwarfs due to the 
water absorption from 2.07 /im to 2.38 //m. 
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Fig. 30. — Comparison of BT-Settl-2010 synthetic spectra by 
Spec K band spectra of HIP 12961, Gl 324 B and LHS 2090. The most prominent molecular 
and atomic features are indicated. The overall shape of the observed spectra, dominated 
mostly by water absorption, are well matched by the BT-Settl-2010 models. The strong 
Na I doublet and the Ca I triplet features in these stars cannot be reproduced by any of the 
supersolar metallicity models. Perhaps the best match is for HIP 12961, where the [M/H] 
=+0.3 and T e g= 4100K BT-Settl-2010 model matches the strengths of the CO bands and 
the neutral Ti and Al lines quite well. The BT-Settl-2010 model modestly underpredicts 
the strengths of the 2.205 Na I and 2.26 Ca I lines, however, and also shows stronger Ca I 
features at ~1.97 /im, and no Mg I line at ~2. 11 fim. Uncertain oscillator stre ngths and 



missing opacity sources could explain these discrepancies ( jRajpurohit et al. 



20101 ). 
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Fig. 31. — Equivalent widths of the Na I doublet (top) and Ca I triplet (middle), and 
H9Q-K2 index (bo ttom), measured from the BT-Settl-2010 synthetic spectra computed by 



Allard et al. 



( 120101 ) and shown as a function of model T e ff- The strengths of the Na I and 
Ca I features are somewhat weaker in the synthetic spectra than in observed data, but the 
qualitative behavior of the Ca I feature as a function of both temperature and metallicity 
is consistent. Water absorption appears to be a monotonic function of temperature, inde- 
pendent of metallicity, for models with T e g- > 3000K. Below 3000K, however, the H2O-K2 
indices measured from the [M/H] < -0.5 dex spectral grids diverge from the [M/H] > -0.5 
models; the stark divergence between the two domains suggest the difference may be a com- 
putational artifact rather than a true astrophysical difference. The Na I doublet shows the 
oddest behavior in these spectra, with significant structure as a function of both metallicity 



and temperature. 
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Fig. 32. — The H 2 0-K2 index measured from the BT-Settl-2010 synthetic spectra shown 
as a function of model T e ff- Triangles represent models with log g = 4.5 and black dots 
represent models with log g= 5.0. The H 2 0-K2 index shows negligible sensitivity to surface 
gravity in all models with T e //> 3000 K. Differences in H 2 0-K2 index due to surface gravity 
are small in the solar and super-solar models throughout the whole effective temperature 
range. The largest discrepancies due to surface gravity are found at lower temperatures in 
the subsolar [M/H] models. The largest discrepancy between the two surface gravities in the 
[M/H]=0.0 dex models is believed to be a computational artifact. 
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Fig. 33. — Top: The H2O-K2 index versus the Spectral Type of M dwarfs. The red dots 
represent the fifty four stars with optical KHM spectral types by RECONS. The dashed 
black line is the linear relationship between KHM spectral types and the H2O-K2 index in 
Equation El obtained using a bayesian approach that includes the errors in both coordinates 
of the objects. Bottom: Spectral type residuals vs H 2 0-K2 index. 
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Fig. 34. — Mean TripleSpec K band spectral type templates. The stars in our sample were 
organized by spectral type and a template was constructed by combining the spectra in each 
subtype. All of the subtypes are constituted by at least three stars, and up to as many as 
thirty, with the exception of the M9 subtype, for which LHS 2924 is the only prototype. The 
most prominent molecular and atomic features are indicated. 
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Fig. 35. — [Fe/H] residuals versus the independent variables of the [Fe/H] fit in Equation [7] 
(Na I /H 2 0-K2, left; Ca I /H 2 0-K2, middle), and as a function of the predicted metallicity 
(right). 
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Fig. 36.— The EW of the Ca I vs the EW of the Na I , both weighted by the H 2 0-K2 
index, for the M dwarfs in our sample. The red and blue dots are the 18 M dwarfs in the 
metallicity calibration sample with [Fe/H]>0.0 and [Fe/H]<0.0, respectively. The black stars 
represent the M dwarf planet hosts, and the yellow dots the rest of the M dwarfs studied in 
this work. The dashed line is an iso-metallicity contour for [Fe/H] =0.0. The dotted lines 
are iso-metallicity contours for [Fe/H] values +0.4, +0.2, -0.2, -0.4, and -0.6 dex from the 
top right corner to the bottom left corner, respectively. The isometallicity contours were 
calculated from Equation [71 In this plane, there is a clear distinction between the M dwarfs 
in the calibration sample with metal-rich and metal-poor FGK-dwarf companions. The M 
dwarf planet hosts all have [Fe/H]>0.0 dex, except Gl 581 and Gl 649, whose [Fe/H] are 
equal to -0.10 dex and -0.04 dex, respectively. The most iron-rich star in the sample is the 
planet host HIP 79431 (+0.46 dex) and the least iron-rich is the flare star V* V1513 Cyg 
(-0.64 dex). 
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Fig. 37. — [[M/H] residuals versus the independent variables of the [M/H] fit in Equation M 
(Na I /H 2 0-K2, left; Ca I /H 2 0-K2, middle), and as a function of the predicted metallicity 
(right). 
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Fig. 38.— The EW of the Ca I vs the EW of the Na I both weighted by the H 2 0-K2 index, 
for the M dwarfs in our sample. The nomenclature is the same as Figure ESJ The dashed 
line is a iso-metallicity contour for [M/H] =0.0. The dotted lines are iso-metallicity contours 
for [M/H] values +0.4, +0.2, -0.2, and -0.4 dex from the top right corner to the bottom 
left corner, respectively. The isometallicity contours were calculated from Equation The 
M dwarf planet hosts all have [M/H]>0.0 dex, except Gl 581 and G1649, whose [M/H] are 
equal to -0.06 dex and -0.02 dex. The most metal-rich star in the sample HIP 79431 has a 
[M/H]=+0.33 dex and the least metal-rich V* V1513 Cyg has [M/H]=-0.45 dex. 
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Fig. 39.— EW of Ca I vs EW o : 
synthetic models by 



Allard et al. 



Na I , weighted by H 2 0-K2 index, for the BT-Settl-2010 
( l2010l ) with T e ff< 3000K. The dashed lines are isometallicity 
contours for [M/H] values, from the top right corner to the bottom left corner, of -0.2, -0.4, 
and -0.6. The isometallicity contours for [M/H] were calculated from Equation [HJ There is a 
clear distinction between the different metallicity models in this plane. The BT-Settl-2010 
synthetic models follow the same direction as the isometallicity contours but their weaker 
strengths of the Na I doublet and Ca I triplet do not allow a direct metallicity comparison 
with empirical data. Subsolar and solar models of temperatures between 3700K and 3000K 
appear as parallel lines in this plane, but the supersolar models intersect the solar model 
at 3000K, due to the odd behavior of the Na I doublet features of the supersolar models 



mentioned in the text 
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Fig. 40. — Differences between the [Fe/H] values predicted for M dwarf stars in our sample 
from the K band [Fe/H] fit in Equation [7] and al) the [Fe/H] from B05, bl) the [Fe/H] from 
JA09, cl) the [Fe/H] from SL10, dl) the [Fe/H] from W09. The B05 and W09 calibrations 
underestimate [Fe/H] values of the supersolar K band [Fe/H] stars. The JA09 calibration 
overestimates the [Fe/H] values of the subsolar K band [Fe/H] stars. The scatter of the 
residuals increases for the subsolar K band [Fe/H] stars when their [Fe/H] are compared 
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Fig. 41. — Difference between the spectral types derived using the H 2 0-K2 index and the 
Ti05 (left), CaH2 (middle), and CaH3 (right) indices versus [M/H]. The Ti05 index predicts 
later subtypes than the H 2 0-K2 index for the solar and metal-rich stars, and earlier subtypes 
to the metal-poor stars. The CaH2 and CaH3 indices show less metallicity sensitivity. The 
CaH2 index predicts slightly later subtypes than the H 2 0-K2 index, while the CaH3 index 
shows an opposite offset, but no bias with metallicity is visible for any of these two CaH 
indices. 
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Fig. 42.— The PMSU Ti02, Ti03, and CaHl indices versus the H 2 0-K2 index . M dwarfs 
with [M/H] > 0.0 dex and [M/H]< 0.0 dex are color-coded with red and blue, respectively. 
The TiO indices correlate with the H 2 0-K2 index, with larger scatter towards early type 
M dwarfs. The CaHl index seems to be more sensitive to metallicity than any of the other 
indices. 
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Fig. 43.— EW of Ca I vs EW of Na I , weighted by H 2 0-K2 index, for the M-M binaries in 
Table [SJ The binary components are linked with lines. The dotted lines are isometallicity 
contours for [M/H] calculated from Equation El The red dots/lines indicate the systems 
where the K band metallicity estimates for the individual components agree to within 0.02 
dex. 
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Fig. 44. — Toomre diagram for the stars in our sample with PMSU kinematics, color-coded 
by metallicity. The metal-poor stars ([Fe/H]<-0.14 dex), depicted by blue dots, cover the 
whole velocity space, while the metal-rich stars ([Fe/H]>0.0 dex) exhibit V to t velocities no 
larger than ~100 km/s. 
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Fig. 45. — K band metallicity distributions for the stars in our sample with PMSU kinemat- 
ics. The boundaries of these kinematic sub-groups are: thin disk, V( O t<50 km/s; thin-thick, 
50 km/s <V toi <70 km/s; thick disk, 70 km/s <V tot <200 km/s. The metallicity distributions 
of these kinematically selected subgroups are consistent with the trends seen for solar-type 
stars: stars with large space motions are preferentially metal-poor, while stars with small 
space motions are preferentially metal-rich. 
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Fig. 46. — H Q st rength as function of spectral type for the stars in our sample with H a 



measurements by 



Gizis et al. 



(120021 ) . The dotted-line represents a 2 A lower limit in H Q 



to define M dwarfs with significantl y enhanced activi ty (i.e., more active than required to 



satisfy the 0.75 A limit adopted by 



West et al. 



(I2008h ). All late-type (>M5) stars are still 



active in all three metallicity categories, but activity persists to earlier types in the metal- 
rich samples, consistent with an underlying correlation between a star's age, activity, and 
metallicity. Chromospheric activity is noticeable for M3 stars with [M/H]>-0.01 dex, for M4 
stars with -0.15<[M/H]<-0.01 dex, and for a M5 stars with [M/H]<-0.15 dex. V* V1513 
Cyg is a metal-poor M2 dwarf that still exhibits chromospheric activity. 
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Fig. 47. — [Fe/H] distribution of 35 M dwarfs whose RV measurements have rule out the 
presence of Jupiter-size planets within several AU by the CPS team, along with the [Fe/H] 
distribution of Jupiter hosts and Neptune hosts in our sample. 
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Fig. 48.— Optical PMSU spectra and K band spectra of LHS 3409 (top), G 203-47 (middle), 
and Gl 169. 1A (bottom). The PMSU TiO and CaH indices are shown in red and blue, 
respectively. The most prominent atomic features in the M dwarf K band spectra are also 
indicated. The regions used to calculate the water index , H 2 0-K2, are shown in gray. G 
203-47 and Gl 169.1 A belong to binary systems with white-dwarf companions. LHS 3409, 
a metal-poor star ([M/H]=-0.36) of similar spectral type, is shown for comparison. The K 
band absorption features are weaker in the spectrum of LHS 3409, but the CaH absorption 
in its optical spectrum is stronger than in the 2 metal-rich stars, which indicates that G 
203-47 and Gl 169. 1A have supersolar metallicities. 



Table 1. TripleSpec Nearby M dwarf Sample - Source Properties 



Name d (pc) V (mag) V Ref. 2MASS K Date SNR Star 



LHS 3576 


19.9 


10.29 




6.553 


2009-07-20 


278 


LSPM 


Gl 338 A 


6.2 


7.64 




3.99 


2009-02-05 


243 


8pc 


Gl 338 B 


6.3 


7.7 




4.14 


2009-02-05 


350 


8pc 


G 210-45 


23.4 


11.21 




7.251 


2009-07-20 


290 


LSPM 


Gl 205 


5.7 


7.97 




4.039 


2009-09-30 


652 


8pc 


Gl 725 A 


3.6 


8.94 




4.432 


2009-07-19 


866 


8pc 


Gl 412 A 


4.8 


8.82 




4.769 


2009-02-05 


337 


8pc 


Gl 686 


8.1 


9.62 




5.572 


2009-07-19 


702 


LSPM 


Gl 752 AB 


5.9 


9.12 




4.673 


2009-07-19 


416 


8pc 


Gl 649 


10.3 


9.7 




5.624 


2010-05-22 


862 


Dlanet/LSPM 


HIP 12961 


23.8 


10.25 




6.736 


2010-11-25 


268 


Dlanet/LSPM 


\ji ziz 


in c 
IZ.O 


y. i o 




o. I oy 


onnn no qh 

zuuy-uy-ou 


OZ / 


caiiDrator/ Ldt ivi 


HD 46375 B 


26.4 


11.8 




7.843 


2009-02-04 


795 


calibrator/LSPM 


HIP 79431 


14.9 


11.34 




6.589 


2010-05-22 


659 


planet/LSPM 


V* V1513 Cyg 




11.98 




8.113 


2009-07-20 


864 


LSPM 


Gl 411 


2.5 


7.49 




3.254 


2009-02-05 


444 


8pc 


Gl 526 


5.4 


8.46 




4.415 


2010-05-22 


1996 


8pc 


V* V547 Cas 


10.1 


10.27 




6.037 


2007-11-17 


385 


LSPM 


Gl 872 B 


16.3 


11.7 


3 


7.3 


2009-07-20 


356 


calibrator/LSPM 



Table 1 — Continued 



Name 


d (pc) 


V (mag) 


V Ref. 


2MASS K 


Date 


SNR 


Star 


Gl 797 B 


20.1 


11.88 


8 


7.416 


2009-07-19 


581 


calibrator /LSPM 


LHS 3577 


12.5 


10.79 


1 


6.533 


2009-07-20 


311 


LSPM 


Gl 581 


6.3 


10.57 


1 


5.837 


2009-02-05 


359 


planet /8pc 


Gl 408 


6.6 


10.03 


1 


5.503 


2009-12-24 


2079 


8pc 


Gl 251 


5.5 


9.89 


1 


5.275 


2009-02-04 


422 


8pc 


Gl 297 2 R 


12.5 


11.8 


3 


7.418 


2008-02-16 


260 


LSPM 


Gl 250 R 


8.7 


10.05 


3 


5.723 


2009-02-04 


681 


calibrator /LSPM 


Gl 176 


9.4 


9.95 


1 


5.607 


2009-09-30 


667 


nlanet /LSPM 

JJldllHj / J_J vJ J- lvl_ 


NLTT 14186 


34.5 


14.60 


8 


7.621 


2009-09-30 


630 


calibrator /LSPM 


Gl 849 


8.8 


10.41 


1 


5.594 


2009-07-20 


755 


nlanet /LSPM 

JJldllHj / J_JkJ-L lv± 


Gl 725 R 


3.5 


9.7 


1 


5 


2009-07-19 


832 


8pc 


Gl 661 AR 


6.1 


9.44 


3 


4.83 


2009-07-19 


708 


8pc 


G 262-29 


32.8 


11.7 


2 


7.61 


2009-07-20 


425 


LSPM 


LHS 3605 


13.7 


11.98 


1 


7.64 


2009-07-20 


307 


LSPM 


LHS 115 


10.2 


12.19 


3 


6.377 


2007-11-17 


326 


LSPM 


Gl 625 


6.6 


10.13 


1 


5.833 


2009-07-19 


573 


8pc 


Gl 643 


6.5 


11.73 


1 


5.756 


2009-07-19 


530 


8pc 


Gl 273 


3.8 


9.84 


1 


4.857 


2009-02-04 


327 


8pc 


LHS 3591 


32.4 


12.73 


2 


8.238 


2009-07-20 


288 


LSPM 



Table 1 — Continued 



Name 


i / \ 

d (pc) 


V (mag) 


V Ref. 


2MASS K 


Date 


SNR 


Star 


Gl 860 AB 


4 


9.59 


1 


4.777 


2009-07-20 


304 


8pc 


Gl 687 


4.5 


9.15 


1 


4.548 


2010-05-24 


984 


8pc 


Gl 628 


4.3 


10.1 


1 


5.075 


2010-05-24 

Z_ V 7 1 V 7 V 7 \ / 1 


808 


8pc 


Gl 873 


5 


10.29 


1 


5.299 


2009-07-20 


392 


8pc 


LHS 3558 


8 


10.54 


1 


5.933 


2009-07-20 


628 


8pc 


G 168-24 


16.3 


12.51 


2 


7.873 


2009-07-20 


752 


LSPM 


HD 222582 R 


41.9 


14.50 


7 


9.583 


2009-07-20 


376 




Gl 436 


10.2 


10.67 


1 


6.073 


2008-02-1 7 


997 


nlanet /LSPM 


HTP 57050 


11 


11.86 


1 


6.822 


2010-05-22 


537 


nlanot /LSPM 

JJldllHj/ 1 J ..J 1 ±V_L 


LP 816-60 

1 J 1 V../ 1 V / V / V 7 


5.5 


11.41 


1 


6.199 


2009-07-19 


899 


8pc 


Gl 896 A 


6.2 


10.05 


1 


5.326 


2009-07-20 


923 


8pc 


Gl 876 


4.7 


10.16 


1 


5.01 


2009-07-20 


577 


TtlflTlpt / ^TiC 


Gl 402 


5.6 


11.64 


1 


6.371 


2009-12-24 


1452 


8pc 


Gl 53,1 B 


20.7 


13.6 


3 


8.673 


2007-11-16 


508 


LSPM 


Gl 555 


6.1 


11.32 


1 


5.939 


2010-05-22 


716 


8pc 


Gl 179 


12.1 


11.94 


1 


6.942 


2010-11-25 


278 


planet/LSPM 


LHS 494 


15.9 


12.51 


2 


7.397 


2009-07-20 


491 


LSPM 


Gl 388 


4.9 


9.4 


2 


4.593 


2009-02-05 


260 


8pc 


Gl 169.1 A 


5.5 


11.06 


2 


5.717 


2009-09-30 


1508 


8pc 



Table 1 — Continued 



Name 


d (pc) 


V (mag) 


V Ref. 


2MASS K 


Date 


SNR 


Star 


LHS 3409 


20.3 


15.11 


2 


10.229 


2010-05-24 


280 


LSPM 


Gl 699 


1.8 


9.54 


I 


4.524 


2010-05-24 


758 


8pc 


LHS 220 


13.3 


13.77 


2 


8.773 


2007-1 1-1 7 


302 


LSPM 


Gl 783 2 B 


20.4 


13.94 


3 


8.883 


2009-07-19 


602 


calibrator /LSPM 


Gl 445 


5.4 


10.8 


1 


5.954 


2009-02-05 


201 


8pc 


Gl 213 


5.8 


11.56 


1 


6.389 


2009-02-04 


318 


8pc 


Gl 544 B 


19 


14.5 


3 


9.592 


2008-02-1 6 


811 


calibrator /LSPM 


LHS 1723 


7.5 


12.16 


3 


6.736 


2009-09-30 


582 


8pc 


GT 1224 


7.5 


13.64 


2 


7.827 


2009-07-19 


876 


8pc 


GJ 1119 


10.3 


13.32 


2 


7.741 


2007-1 1-1 6 


972 


LSPM 


G 041-014 


5.6 


10.89 


3 


5.688 


2009-12-24 


857 


8pc 


GT 3348 B 


23.5 


13.98 


3 


8.791 


2009-09-30 


416 


calibrator /LSPM 


t lit' -i o/~\r\ 

LHb 1809 


9.3 


14.45 


2 


8.435 


2007-11-17 


746 


LSPM 


Gl 447 


3.3 


11.12 


1 


5.654 


2009-02-05 


304 


8pc 


LHS 1066 


16.6 


14 


3 


9.11 


2007-11-17 


668 


LSPM 


GJ 3134 


24 


14.31 


3 


8.995 


2007-11-17 


499 


LSPM 


Gl 231.1 B 


19.9 


13.27 


3 


8.267 


2009-02-04 


908 


calibrator/LSPM 


LHS 3593 


13.9 


13.98 


2 


8.481 


2009-07-20 


228 


LSPM 


GJ 3379 


5.4 


11.33 


2 


6.042 


2009-02-04 


662 


8pc 



Table 1 — Continued 



Name 


d (pc) 


V (mag) 


V Ref. 


2MASS K 


Date 


SNR 


Star 


Gl 268 AB 


6.4 


11.65 


1 


5.846 


2009-02-04 


443 


8pc 


ni 768 1 R 


1 9 4 


1 3 1 


3 


8 01 2 


201 0-05-22 


531 

Jul 


calibrator /LSPM 


NLTT 25869 


11 


14.5 


3 


8.64 


2008-02-16 


551 


LSPM 


LHS 6007 


21.3 


14.25 


2 


8.852 


2007-11-17 


437 


LSPM 


Gl 905 


3.2 


12.27 


2 


5.929 


2009-07-20 


824 


8pc 


LHS 495 


Q 8 


13.41 


2 


7 749 


2009-07-20 


258 


LSPM 


GJ 1214 


13 


14.67 


2 


8.78 


2010-05-24 


415 


olanet/LSPM 


G 246-33 


14 


14.63 


3 


8.656 


2007-1 1-1 6 


618 


LSPM 


Gl 324 R 


12.5 


13.15 


3 


7.666 


2008-02-1 6 


814 


calibrator /LSPM 

LoJlUI OitUl / 1 j kJ 1 ±V_L 


G 203-47 


7.3 


11.77 


1 


6.485 


2009-07-1 9 


659 


8pc 


Gl 299 


6.8 


12.83 


3 


7.66 


2009-02-04 


858 


8pc 


LHS 224 


9.2 


13.3 


2 


7.776 


2007-1 1-1 7 


465 


LSPM 


Gl 611 R 


13.8 


14.2 


2 


9.159 


2008-02-16 


970 


calibrator/LSPM 


Gl 630.1 A 


14.5 


12.9 


2 


7.796 


2008-02-17 


483 


LSPM 


NSV 13261 


15.9 


14.67 


2 


8.753 


2009-07-20 


407 


LSPM 


NLTT 15867 


25 


16.49 


3 


10.312 


2009-02-04 


334 


calibrator/LSPM 


Gl 166 C 


5 


11.17 


3 


5.962 


2009-09-30 


591 


calibrator/8pc 


LHS 3376 


7.3 


13.46 


2 


7.948 


2009-07-19 


718 


8pc 


GJ 3253 


14.3 


13.79 


2 


8.053 


2007-11-16 


417 


LSPM 



Table 1 — Continued 



Name 


d (pc) 


V (mag) 


V Ref. 


2MASS K 


Date 


SNR 


Star 


GT 3069 


14.4 


15.12 


3 


8.864 


2007-11-17 

iiUU 1 J- J- -LI 


412 


LSPM 


Gl 1986 


7 9 


1 4 69 


2 


8 183 

O. -LOO 


9009-07-90 

ZjVJVJJ? VJ 1 Zj\J 


228 




Gl 164 

V. 1 1 1_ VJ4: 


1 1 Q 


1 3 5 

lO . O 


2 


7 915 


9008-09-1 7 

ZjVJVJO VJZj -L 1 


495 


LSPM 


Gl 777 B 


11.7 


14.33 


2 


8.712 


9009-07-19 

ZiVJVJi? VJ 1 li7 


522 


ralihrator /LSPM 


Gl 866 


3.4 


12.18 


3 


5.537 


9009-07-90 

ZiVJVJi? VJ 1 Zj VJ 


400 


8pc 


Gl 473 AR 


4.3 


12.44 


3 


6.042 


9009-09-05 

ZiVjvji? vj Zj vjfj 


413 


8pc 


Gl 934 AR 


4.1 


11.12 


I 


5 486 


9009-09-04 

ZjUVJC/ VJZj VJt: 


691 




LHS 3549 


9.3 


14.04 


2 


8.095 


9009-07-90 

Z/VJVJi7 VJ I Zj VJ 


442 


LSPM 


LHS 1706 


14.1 


15.23 


2 


8.977 


2009-09-30 


372 


LSPM 


V* V388 Gas 

V V OOO V_ ClO 




13.78 


5 


7.718 


9007-1 1-1 7 

ZjVJVJ I _L _L _L ( 


410 


LSPM 


LHS 37QQ 


7.4 


13.25 


2 


7.319 


9009-07-90 

ZiVJVJC/ VJ I Zj VJ 


388 


8pc 


Gl 285 


5.9 


11.19 


1 


5.698 


9009-09-04 


1447 


8pc 


LHS 18 


12 


14.66 


2 


9.083 


2007-11-17 


475 


LSPM 


Gl 412 R 


4.8 


14.45 


3 


7.839 


2009-02-05 


367 


8pc 


LHS 1901 


8 


15.87 


3 


9.126 


2007-11-17 


538 


8pc 


LSPM J0011 a 


11.7 


15.87 


3 


9.093 


2007-11-17 


666 


LSPM 


GJ 1245 AC 


4.5 


13.41 


2 


6.854 


2009-07-19 


539 


8pc 


GJ 1245 R 


4.7 


13.99 


3 


7.387 


2009-07-19 


440 


8pc 


GJ 1116 AR 


5.2 


14.06 


2 


6.889 


2009-02-05 


1094 


8pc 



Table 1 — Continued 



Name 


d (pc) 


it/ \ 

V (mag) 


V Rcf. 


2MASS K 


Date 


SNR 


Star 


GT 3146 


8.5 


15.79 


2 


8.981 


2007-11-17 


485 


LSPM 


LHS 252 


10 


15.05 


2 


8.668 


2007-1 1-1 7 


772 


LSPM 


C\ 376 B 


1 3 Q 


16 13 

1U . J- O 


Q 


9 275 


2009-1 2-24 


372 


ralihratnr /LSPM 


Gl 1156 


6.5 


13.79 


2 


7.57 


2009-02-05 


505 


8pc 


Gl 406 


2.4 


13.54 




6.084 


2009-1 2-24 


1860 


8pc 


GJ 3147 


10.4 


15.99 


2 


9.011 


2007-11-16 


507 


LSPM 


LHS 2Q2 


4 5 


1 5 6 


2 


7 926 


2009-1 2-24 


870 




Gl 644 C 


6.5 


16.7 




8.816 


2009-07-19 


542 


8pc 


GJ 1111 


3.6 


14.81 


2 


7.26 


2008-02-16 


746 


8pc 


LHS 20Q0 


q 


16.1 


3 


8.437 


2009-1 2-24 


914 


8pc 


V* V4Q2 Lvr 

V V J_i V ± 


14.1 


18.23 


2 


10.308 


2009-07-20 


370 


LSPM 


'Tipper n r H pn 


3.6 


15.4 


3 


7.585 


2009-02-04 


1092 


8pc 


2MASS J18 


6.2 


18.27 


3 


9.171 


2009-07-19 


383 


8pc 


LHS 2065 


8.5 


18.8 


2 


9.942 


2008-02-16 


264 


LSPM 


LHS 2924 


10.8 


19.74 


2 


10.744 


2008-02-17 


560 


LSPM 


Gl 908 c 


6.0 


8.98 


3 


5.043 


2009-07-20 


253 


8pc 


Gl 809 c 


7.0 


8.54 


3 


4.618 


2009-07-19 


680 


8pc 


Gl 644 AB c 


6.5 


9.80 


3 


4.403 


2009-07-19 


587 


8pc 


Gl 829 AB c 


6.7 


10.30 


3 


5.453 


2009-07-19 


208 


8pc 



Table 1 — Continued 



Name d (pc) V (mag) V Ref. 2MASS K Date SNR Star 



a LSPM J0011+5908 

b 2MASS J18353790+3259545 

c Stars with low quality K band spectra 



*Reference of di stance and V magnitu de: (1) HIPP ARCOS (Ivan Leeuwen 



20071). (2) YALE 



Koenetal 



20041) . (8) 



van Altena et al 



(201 01), (5) 



Gliese fc Jahreissl ( 199ll) 



Leggetd (119921) . (6) 



2001 ), (3) PMSU dReid et al 



Bessell(ll990l) . (7) 



1997 ), (4) 



Gould fc Chaname 
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Table 2. EW Bandpasses and Continuum Points 



Feature 


Wavelength 


Integration Limits 


Continuum Points 




H 


H 


H 


Na I 


2.206 


2.2020-2.2120 


2.1965, 2.2125,2.2175 




2.209 






Ca I 


2.261 


2.2580-2.2690 


2.2510, 2.2580, 2.2705, 2.2760 



2.263 
2.265 



Table 3. TripleSpec Nearby M dwarf Sample - Spectral Measurements 



Name EW(Na I) EW(Ca I) H 2 0-K2 Sp.T. T eff [M/H] [Fc/H] 



LHS 3576 


4 


.355 


± 





.107 


3 


.883 


± 





095 


1 


.021 


± 





.004 





3883 ± 24 


-0.13 


± 


0.12 


-0.19 


± 


0.17 


Gl 338 A 


4 


.045 


± 





.082 


4. 


.358 


± 


0. 


.081 


1 


.038 


± 





.003 





4031 ± 56 


-0.13 


± 


0.12 


-0.18 


± 


0.17 


Gl 338 B 


4 


.023 


± 





.119 


4. 


475 


± 





.107 


1 


.019 


± 





.004 





3869 ± 15 


-0.11 


± 


0.12 


-0.15 


± 


0.17 


G 210-45 


4 


.998 


± 





104 


4. 


.633 


± 


0. 


097 


1 


.024 


± 





.004 





3910 ± 37 


-0.03 


± 


0.12 


-0.05 


± 


0.17 


Gl 205 


7 


.636 


± 





.043 


6 


105 


± 


0. 


043 


1 


.032 


± 


0. 


002 





4012 ± 106 


0.25 


± 


0.12 


0.35 


± 


0.17 


Gl 725 A 


3 


.068 


± 





.034 


2 


152 


± 


0. 


031 





.982 


± 





.002 


1 


3680 


± 


18 


-0.34 


± 


0.12 


-0.49 


± 


0.17 


Gl 412 A 


3 


.431 


± 





091 


2 


.618 


± 


0. 


.085 





.983 


± 





.003 


1 


3684 


± 


20 


-0.28 


± 


0.12 


-0.40 


± 


0.17 


Gl 686 


3 


.855 


± 





039 


3 


.268 


± 


0. 


037 





.985 


± 


0. 


002 


1 


3693 


± 


20 


-0.20 


± 


0.12 


-0.28 


± 


0.17 


Gl 752 AB 


5 


.626 


± 





074 


4. 


.018 


± 


0. 


068 


1 


.003 


± 


0. 


.003 


1 


3789 


± 


20 


-0.03 


± 


0.12 


-0.05 


± 


0.17 


Gl 649 


4 


.966 


± 


0. 


.087 


4. 


534 


± 


0. 


033 





.993 


± 


0. 


002 


1 


3733 


± 


20 


-0.02 


± 


0.12 


-0.04 


± 


0.17 


HIP 12961 


5 


.175 


± 





105 


4. 


963 


± 


0. 


105 


1 


.014 


± 





004 


1 


3838 


± 


19 


0.01 


± 


0.12 


0.01 


± 


0.17 


Gl 212 


5 


.078 


± 





045 


5. 


.206 


± 


0. 


045 


1 


.016 


± 


0. 


002 


1 


3851 


± 


17 


0.02 


± 


0.12 


0.03 


± 


0.17 


HD 46375 B 


7 


.136 


± 


0. 


035 


5 


.407 


± 





034 





.978 


± 


0. 


.002 


1 


3663 


± 


15 


0.21 


± 


0.12 


0.29 


± 


0.17 


HIP 79431 


8 


.948 


± 


0. 


041 


5 


.436 


± 





039 





.984 


± 


0. 


.002 


1 


3689 


± 


20 


0.33 


± 


0.12 


0.46 


± 


0.17 


V* V1513 Cyg 


1 


.508 


± 





.061 


1 


.964 


± 





.059 





.935 


± 





.002 


2 


3483 


± 


17 


-0.45 


± 


0.12 


-0.64 


± 


0.17 


Gl 411 


3 


.115 


± 





065 


2 


.623 


± 


0. 


059 





.945 


± 





.002 


2 


3526 


± 


18 


-0.28 


± 


0.12 


-0.41 


± 


0.17 


Gl 526 


3 


.646 


± 


0. 


015 


3 


207 


± 





014 





.973 


± 


0. 


.001 


2 


3642 


± 


17 


-0.21 


± 


0.12 


-0.30 


± 


0.17 


V* V547 Cas 


3 


.603 


± 





.075 


3 


.246 


± 


0. 


068 





.945 


± 


0. 


003 


2 


3526 


± 


18 


-0.19 


± 


0.12 


-0.28 


± 


0.17 


Gl 872 B 


4 


101 


± 


0. 


.084 


3 


173 


± 


0. 


075 





.955 


± 


0. 


003 


2 


3569 


± 


20 


-0.17 


± 


0.12 


-0.25 


± 


0.17 


Gl 797 B 


4 


.298 


± 





.054 


3 


.182 


± 


0. 


045 





.955 


± 





002 


2 


3569 


± 


20 


-0.16 


± 


0.12 


-0.23 


± 


0.17 


LHS 3577 


4 


.704 


± 





.095 


3 


.919 


± 


0. 


085 





.968 


± 


0. 


003 


2 


3621 


± 


19 


-0.07 


± 


0.12 


-0.11 


± 


0.17 


Gl 581 


4 


.965 


± 





078 


3 


.656 


± 





077 





.947 


± 


0. 


.003 


2 


3534 


± 


18 


-0.06 


± 


0.12 


-0.10 


± 


0.17 


Gl 408 


4 


.755 


± 





.015 


3 


.895 


± 


0. 


013 





.945 


± 


0. 


001 


2 


3526 


± 


18 


-0.06 


± 


0.12 


-0.09 


± 


0.17 


Gl 251 


4 


.935 


± 





.047 


3 


862 


± 


0. 


.053 





.939 


± 





.002 


2 


3500 


± 


17 


-0.04 


± 


0.12 


-0.07 


± 


0.17 



Table 3 — Continued 



Name EW(Nal) EW(Ca I) H 2 0-K2 Sp.T. T eff [M/H] [Fe/H] 



Gl 297.2 B 


5 


.154 


± 





.107 


4. 


.055 


± 





.125 





.959 


± 


0. 


.003 


2 


3585 


± 


20 


-0.03 ± 


0.12 


-0.04 ± 


0.17 


Gl 250 B 


5 


.124 


± 





.039 


4. 


.427 


± 





.041 





.955 


± 





.002 


2 


3569 


± 


20 


0.01 ± 


0.12 


0.01 ± 


0.17 


Gl 176 


6 


.155 


± 


0. 


.042 


4. 


.806 


± 





.042 





.958 


± 





.002 


2 


3581 


± 


20 


0.11 ± 


0.12 


0.15 ± 


0.17 


NLTT 14186 


6 


.259 


± 





047 


5. 


.302 


± 





.041 





.946 


± 





.002 


2 


3530 


± 


18 


0.17 ± 


0.12 


0.24 ± 


0.17 


Gl 849 


7. 


.198 


± 


0. 


055 


5, 


382 


± 





.049 





963 


± 


0. 


002 


2 


3601 


± 


19 


0.23 ± 


0.12 


0.31 ± 


0.17 


Gl 725 B 


3 


594 


± 


0. 


035 


2. 


.308 


± 





.033 





.895 


± 


0. 


002 


3 


3288 


± 


27 


-0.25 


± 


0. 


12 


-0.36 


± 


0.17 


Gl 661 AB 


3 


.577 


± 


0. 


039 


2. 


714 


± 


0. 


035 


0. 


.892 


± 





.002 


3 


3272 


± 


28 


-0.21 


± 


0. 


12 


-0.31 


± 


0.17 


G 262-29 


3 


.474 


± 


0. 


072 


3. 


.056 


± 





.065 





,929 


± 





.002 


3 


3455 


± 


16 


-0.21 


± 


0. 


12 


-0.30 


± 


0.17 


LHS 3605 


3 


.946 


± 


0. 


.098 


2. 


879 


± 





.086 





.925 


± 





.003 


3 


3437 


± 


16 


-0.19 


± 


0. 


12 


-0.28 


± 


0.17 


LHS 115 


4. 


101 


± 


0. 


.085 


2 


.683 


± 





.078 





.905 


± 


0. 


.004 


3 


3339 


± 


23 


-0.19 


± 


0. 


12 


-0.27 


± 


0.17 


Gl 625 


4. 


.112 


± 





.049 


3 


.144 


± 





.046 





.931 


± 





.002 


3 


3464 


± 


17 


-0.16 


± 


0. 


12 


-0.23 


± 


0.17 


Gl 643 


4. 


.703 


± 


0. 


.055 


2 


.617 


± 





.052 





.912 


± 





.002 


3 


3376 


± 


20 


-0.15 


± 


0. 


12 


-0.22 


± 


0.17 


Gl 273 


4. 


.667 


± 


0. 


.085 


2 


915 


± 





.104 





.896 


± 


0. 


.003 


3 


3293 


± 


26 


-0.12 


± 


0. 


12 


-0.17 


± 


0.17 


LHS 3591 


4. 


.777 


± 


0. 


.104 


3. 


.416 


± 





.092 





.923 


± 


0. 


.003 


3 


3428 


± 


16 


-0.08 


± 


0. 


12 


-0.12 


± 


0.17 


Gl 860 AB 


5 


.075 


± 


0. 


097 


3 


.348 


± 





.088 





.931 


± 





.003 


3 


3464 


± 


17 


-0.07 


± 


0. 


12 


-0.11 


± 


0.17 


Gl 687 


4. 


.643 


± 





.057 


3 


.712 


± 





.027 





.916 


± 





.002 


3 


3395 


± 


18 


-0.06 


± 


0. 


12 


-0.09 


± 


0.17 


Gl 628 


5 


.548 


± 


0. 


036 


3 


508 


± 





.034 





.913 


± 





.002 


3 


3380 


± 


20 


-0.01 


± 


0. 


12 


-0.02 


± 


0.17 


Gl 873 


5. 


.573 


± 





.072 


3 


.642 


± 





.068 





.917 


± 


0. 


.003 


3 


3400 


± 


18 


-0.00 


± 


0. 


12 


-0.01 


± 


0.17 


LHS 3558 


4. 


.927 


± 





.048 


4. 


.312 


± 





.042 





.925 


± 


0. 


.002 


3 


3437 


± 


16 


0.01 ± 


0.12 


0.01 ± 


0.17 


G 168-24 


5 


.266 


± 


0. 


.043 


4. 


.029 


± 





.037 





.923 


± 





.002 


3 


3428 


± 


16 


0.01 ± 


0.12 


0.01 ± 


0.17 


HD 222582 B 


5 


.622 


± 





.076 


3 


.607 


± 





.075 





.895 


± 





.003 


3 


3288 


± 


27 


0.02 ± 


0.12 


0.02 ± 


0.17 


Gl 436 


5 


.373 


± 


0. 


.029 


4. 


302 


± 





.024 





.932 


± 





.002 


3 


3469 


± 


17 


0.03 ± 


0.12 


0.04 ± 


0.17 


HIP 57050 


5. 


.659 


± 





.053 


3 


806 


± 





.051 





.895 


± 


0. 


.002 


3 


3288 


± 


27 


0.04 ± 


0.12 


0.05 ± 


0.17 


LP 816-60 


6 


.185 


± 


0. 


034 


3 


.725 


± 





.031 





.918 


± 





.002 


3 


3405 


± 


17 


0.05 ± 


0.12 


0.06 ± 


0.17 



Table 3 — Continued 



Name EW(Nal) EW(Ca I) H 2 0-K2 Sp.T. T e// [M/H] [Fe/H] 



Gl 896 A 


5 


861 


± 





032 


4 


.515 


± 





.026 





.904 


± 





.002 


3 


3334 


± 


23 


0.11 ± 


0.12 


0.15 


± 


0.17 


Gl 876 


6 


777 


± 


0. 


049 


4 


.363 


± 





.047 





.933 


± 





.002 


3 


3473 


± 


17 


0.14 ± 


0.12 


0.19 


± 


0.17 


Gl 402 


7. 


144 


± 


0. 


021 


3 


.882 


± 


0. 


017 


0. 


.904 


± 





.001 


3 


3334 


± 


23 


0.15 ± 


0.12 


0.20 


± 


0.17 


Gl 53.1 B 


7. 


.006 


± 


0. 


055 


3 


.954 


± 





.055 





.894 


± 





.002 


3 


3282 


± 


27 


0.15 ± 


0.12 


0.21 


± 


0.17 


Gl 555 


7. 


.051 


± 


0. 


041 


3 


.992 


± 





.039 





.895 


± 





.002 


3 


3288 


± 


27 


0.16 ± 


0.12 


0.22 


± 


0.17 


Gl 179 


6 


793 


± 


0. 


105 


4 


.554 


± 





.098 


0. 


.922 


± 





.003 


3 


3424 


± 


16 


0.17 ± 


0.12 


0.23 


± 


0.17 


LHS 494 


6 


477 


± 


0. 


061 


4 


.619 


± 


0. 


054 


0. 


.895 


± 





.002 


3 


3288 


± 


27 


0.17 ± 


0.12 


0.24 


± 


0.17 


Gl 388 


6 


.523 


± 


0. 


107 


5 


.105 


± 





.098 





.915 


± 





.003 


3 


3390 


± 


19 


0.20 ± 


0.12 


0.28 


± 


0.17 


Gl 169.1 A 


8. 


323 


± 


0. 


019 


4 


.058 


± 


0. 


018 





.892 


± 





.001 


3 


3272 


± 


28 


0.26 ± 


0.12 


0.36 


± 


0.17 


LHS 3409 


3 


254 


± 


0. 


.235 


1 


.313 


± 





.109 





.875 


± 





.003 


4 


3178 


± 


41 


-0.36 


± 


0. 


12 


-0.52 


± 


0.18 


Gl 699 


3 


745 


± 


0. 


043 


1 


.967 


± 


0. 


037 


0. 


891 


± 





.002 


4 


3266 


± 


29 


-0.27 


± 


0. 


12 


-0.39 


± 


0.17 


LHS 220 


3 


.905 


± 


0. 


094 


1 


.966 


± 





.091 





.863 


± 





.003 


4 


3093 


± 


70 


-0.24 


± 


0. 


12 


-0.35 


± 


0.17 


Gl 783.2 B 


4. 


032 


± 


0. 


053 


2 


.365 


± 


0. 


044 


0. 


875 


± 





.002 


4 


3178 


± 


41 


-0.20 


± 


0. 


12 


-0.29 


± 


0.17 


Gl 445 


3 


.884 


± 


0. 


145 


2 


.825 


± 





135 





878 


± 





.004 


4 


3195 


± 


38 


-0.17 


± 


0. 


12 


-0.25 


± 


0.17 


Gl 213 


4. 


.358 


± 


0. 


088 


2 


.443 


± 


0. 


.086 





874 


± 





.003 


4 


3167 


± 


47 


-0.17 


± 


0. 


12 


-0.25 


± 


0.17 


Gl 544 B 


5 


263 


± 


0. 


037 


2 


.748 


± 


0. 


.038 





.851 


± 





.002 


4 


3058 


± 


65 


-0.05 


± 


0. 


12 


-0.09 


± 


0.17 


LHS 1723 


5 


.463 


± 


0. 


051 


2 


.837 


± 





.047 





.855 


± 





.002 


4 


3054 


± 


69 


-0.03 


± 


0. 


12 


-0.06 


± 


0.17 


GJ 1224 


6 


032 


± 





.035 


2 


.702 


± 





.033 





.887 


± 





.002 


4 


3245 


± 


31 


-0.03 


± 


0. 


12 


-0.05 


± 


0.17 


GJ 1119 


5. 


,599 


± 





.032 


3 


.091 


± 





.028 





.882 


± 





.002 


4 


3218 


± 


34 


-0.02 


± 


0. 


12 


-0.04 


± 


0.17 


G 041-014 


5 


.324 


± 





033 


3 


.309 


± 





.031 





.875 


± 





.002 


4 


3178 


± 


41 


-0.02 


± 


0. 


12 


-0.03 


± 


0.17 


GJ 3348 B 


5 


.438 


± 


0. 


.073 


3 


.323 


± 





.064 





.877 


± 





.002 


4 


3189 


± 


39 


-0.01 


± 


0. 


12 


-0.02 


± 


0.17 


LHS 1809 


6 


.476 


± 


0. 


.037 


2 


.359 


± 





.036 





.855 


± 





.002 


4 


3054 


± 


69 


-0.00 


± 


0. 


12 


-0.01 


± 


0.17 


Gl 447 


5. 


.655 


± 





.092 


3 


.006 


± 





.088 





.853 


± 





.003 


4 


3065 


± 


69 


-0.00 


± 


0. 


12 


-0.01 


± 


0.17 


LHS 1066 


5 


.459 


± 


0. 


049 


3 


.244 


± 





.048 





.851 


± 





.002 


4 


3058 


± 


65 


0.01 ± 


0.12 


0.00 


± 


0.17 



Table 3 — Continued 



Name EW(Nal) EW(Ca I) H 2 0-K2 Sp.T. T e// [M/H] [Fc/H] 



GJ 3134 


5 


.562 


± 


0. 


.058 


3 


.274 


± 





.055 





.861 


± 





.002 


4 


3080 


± 


74 


0.01 


± 


0.12 


0.01 


± 


0.17 


Gl 231.1 B 


5, 


.602 


± 


0. 


.034 


3 


.458 


± 





.032 





877 


± 





.002 


4 


3189 


± 


39 


0.02 


± 


0.12 


0.02 


± 


0.17 


LHS 3593 


6. 


126 


± 





125 


3 


.095 


± 





123 





868 


± 





.004 


4 


3129 


± 


58 


0.03 


± 


0.12 


0.03 


± 


0.17 


GJ 3379 


6. 


106 


± 


0. 


042 


3 


.502 


± 


0. 


045 


0. 


875 


± 





.002 


4 


3178 


± 


41 


0.06 


± 


0.12 


0.08 


± 


0.17 


Gl 268 AB 


6. 


753 


± 


0. 


.047 


3 


175 


± 


0. 


054 


0. 


876 


± 





.002 


4 


3184 


± 


40 


0.08 


± 


0.12 


0.10 


± 


0.17 


Gl 768.1 B 


5. 


847 


± 


0. 


.056 


4. 


.005 


± 


0. 


051 


0. 


887 


± 





.002 


4 


3245 


± 


.31 


0.08 


± 


0.12 


0.10 


± 


0.17 


NLTT 25869 


6. 


452 


± 


0. 


051 


3 


.229 


± 


0. 


049 


0. 


853 


± 


0. 


002 


4 


3065 


± 


69 


0.08 


± 


0.12 


0.11 


± 


0.17 


LHS 6007 


6. 


546 


± 


0. 


062 


3 


951 


± 


0. 


062 


0. 


.891 


± 





.002 


4 


3266 


± 


29 


0.12 


± 


0.12 


0.16 


± 


0.17 


Gl 905 


7. 


796 


± 


0. 


035 


2 


.732 


± 


0. 


035 


0. 


851 


± 


0. 


002 


4 


3058 


± 


65 


0.14 


± 


0.12 


0.19 


± 


0.17 


LHS 495 


7. 


259 


± 





115 


3 


354 


± 


0. 


104 


0. 


861 


± 


0. 


.003 


4 


3080 


± 


74 


0.15 


± 


0.12 


0.20 


± 


0.17 


GJ 1214 


7. 


753 


± 


0. 


096 


3 


275 


± 


0. 


067 


0. 


887 


± 


0. 


003 


4 


3245 


± 


31 


0.15 


± 


0.12 


0.20 


± 


0.17 


G 246-33 


6. 


468 


± 


0. 


051 


4. 


.232 


± 


0. 


052 





861 


± 





.002 


4 


3080 


± 


74 


0.17 


± 


0.12 


0.24 


± 


0.17 


Gl 324 B 


7. 


818 


± 


0. 


035 


4. 


019 


± 


0. 


034 


0. 


889 


± 


0. 


.002 


4 


3256 


± 


29 


0.22 


± 


0.12 


0.31 


± 


0.17 


G 203-47 


7. 


.978 


± 


0. 


045 


3 


.969 


± 


0. 


043 


0. 


882 


± 


0. 


002 


4 


3218 


± 


34 


0.24 


± 


0.12 


0.33 


± 


0.17 


Gl 299 


3. 


354 


± 


0. 


035 


1 


456 


± 


0. 


032 


0. 


839 


± 


0. 


002 


5 


3021 


± 


49 


-0.32 


± 


0.12 


-0.46 ± 


0.17 


LHS 224 


2. 


905 


± 





063 


1 


762 


± 


0. 


056 


0. 


.826 


± 





.002 


5 


2982 


± 


34 


-0.32 


± 


0.12 


-0.46 ± 


0.17 


Gl 611 B 


3 


277 


± 


0. 


037 


1 


.635 


± 


0. 


029 


0. 


849 


± 


0. 


002 


5 


3051 


± 


61 


-0.32 


± 


0.12 


-0.45 ± 


0.17 


Gl 630.1 A 


3. 


295 


± 


0. 


071 


1 


978 


± 


0. 


067 


0. 


.829 


± 





.002 


5 


2991 


± 


37 


-0.27 


± 


0.12 


-0.39 ± 


0.17 


NSV 13261 


4, 


375 


± 


0. 


075 


1 


707 


± 


0. 


075 


0. 


812 


± 





.003 


5 


2939 


± 


20 


-0.20 


± 


0.12 


-0.29 ± 


0.17 


NLTT 15867 


4 


.513 


± 


0. 


087 


2 


.242 


± 


0. 


.086 


0. 


821 


± 





.003 


5 


2967 


± 


28 


-0.14 


± 


0.12 


-0.21 ± 


0.17 


Gl 166 C 


4. 


522 


± 


0. 


.046 


2. 


.658 


± 


0. 


.049 


0. 


.825 


± 





.002 


5 


2979 


± 


32 


-0.10 


± 


0.12 


-0.15 ± 


0.17 


LHS 3376 


5 


.288 


± 


0. 


.042 


2. 


.333 


± 


0. 


.039 





836 


± 





.002 


5 


3012 


± 


45 


-0.08 


± 


0.12 


-0.12 ± 


0.17 


GJ 3253 


5 


.315 


± 





.068 


2. 


.723 


± 





.065 





.842 


± 





.003 


5 


3030 


± 


52 


-0.05 


± 


0.12 


-0.07 ± 


0.17 


GJ 3069 


5. 


.675 


± 


0. 


.069 


2 


.448 


± 





.065 





839 


± 





.003 


5 


3021 


± 


49 


-0.04 


± 


0.12 


-0.07 ± 


0.17 



Table 3 — Continued 



Name EW(Na I) EW(Ca I) H 2 0-K2 Sp.T. T eff [M/H] [Fe/H] 



GJ 1286 


6 


408 


± 


0. 


.126 


2 


.083 


± 





.129 





.836 


± 


0. 


.004 


5 


3012 


± 


45 


-0.02 


± 


0.12 


-0.04 


± 


0.17 


Gl 164 


5 


.956 


± 


0. 


.055 


2 


.927 


± 





.053 





.849 


± 





.002 


5 


3051 


± 


61 


0.02 


± 


0.12 


0.02 


± 


0.17 


Gl 777 B 


5. 


.857 


± 


0. 


.055 


3 


.005 


± 





.048 





.839 


± 





.002 


5 


3021 


± 


49 


0.03 


± 


0.12 


0.03 


± 


0.17 


Gl 866 


6 


.938 


± 


0. 


.071 


2 


.085 


± 





.065 





.816 


± 


0. 


.002 


5 


2952 


± 


23 


0.04 


± 


0.12 


0.05 


± 


0.17 


Gl 473 AB 


6 


.368 


± 


0. 


.072 


2 


.652 


± 





.085 





.825 


± 


0. 


.003 


5 


2979 


± 


32 


0.05 


± 


0.12 


0.05 


± 


0.17 


Gl 234 AB 


6 


.419 


± 


0. 


.038 


3 


.348 


± 





.039 





.845 


± 


0. 


.002 


5 


3039 


± 


56 


0.10 


± 


0.12 


0.13 


± 


0.17 


LHS 3549 


7. 


.097 


± 





087 


2 


.958 


± 





.059 





.836 


± 





.002 


5 


3012 


± 


45 


0.12 


± 


0.12 


0.16 


± 


0.17 


LHS 1706 


7. 


.681 


± 


0. 


.075 


2 


.822 


± 





.079 





.834 


± 


0. 


.003 


5 


3006 


± 


43 


0.16 


± 


0.12 


0.21 


± 


0.17 


V* V388 Cas 


7. 


.326 


± 


0. 


.067 


3 


.737 


± 





.066 





.845 


± 





.002 


5 


3039 


± 


56 


0.21 


± 


0.12 


0.28 


± 


0.17 


LHS 3799 


8. 


.478 


± 





.075 


3 


.508 


± 





.068 





.845 


± 


0. 


.003 


5 


3039 


± 


56 


0.28 


± 


0.12 


0.38 


± 


0.17 


Gl 285 


7. 


.343 


± 


0. 


.019 


4 


.545 


± 





.019 





.847 


± 





.002 


5 


3045 


± 


59 


0.29 


± 


0.12 


0.40 


± 


0.17 


LHS 18 


3 


285 


± 


0. 


.095 


1 


.213 


± 





.071 





.802 


± 


0. 


.002 


6 


2909 


± 


17 


-0.34 


± 


0.12 


-0.48 


± 


0.17 


Gl 412 B 


5 


.138 


± 





.076 





.385 


± 





.074 





.804 


± 





.003 


6 


2915 


± 


18 


-0.27 


± 


0.12 


-0.39 


± 


0.17 


LHS 1901 


4. 


.322 


± 


0. 


054 


1 


.231 


± 





.053 





.775 


± 





.003 


6 


2826 


± 


21 


-0.23 


± 


0.12 


-0.34 


± 


0.17 


LSPM J0011 a 


5 


.109 


± 





.049 


1 


.285 


± 





.051 





.805 


± 





.003 


6 


2918 


± 


18 


-0.18 


± 


0.12 


-0.27 


± 


0.17 


GJ 1245 AC 


5. 


.595 


± 


0. 


.092 


1 


.675 


± 





.108 





.795 


± 


0. 


.002 


6 


2890 


± 


19 


-0.09 


± 


0.12 


-0.14 


± 


0.17 


GJ 1245 B 


5 


.695 


± 


0. 


068 


1 


.615 


± 





.065 





.792 


± 





.003 


6 


2881 


± 


20 


-0.09 


± 


0.12 


-0.13 


± 


0.17 


GJ 1116 AB 


5 


.735 


± 





.028 


1 


.709 


± 





.025 





.797 


± 


0. 


.002 


6 


2896 


± 


18 


-0.08 


± 


0.12 


-0.12 


± 


0.17 


GJ 3146 


6 


.066 


± 





.062 


1 


.457 


± 





.058 





.788 


± 


0. 


.002 


6 


2868 


± 


21 


-0.07 


± 


0.12 


-0.11 


± 


0.17 


LHS 252 


6 


.004 


± 


0. 


.036 


2 


.287 


± 





.036 





.798 


± 





.002 


6 


2898 


± 


17 


0.00 


± 


0.12 


-0.01 


± 


0.17 


Gl 376 B 


7. 


.725 


± 


0. 


.077 


1 


.687 


± 





.081 





.776 


± 


0. 


.003 


6 


2829 


± 


21 


0.11 


± 


0.12 


0.14 


± 


0.17 


Gl 1156 


6 


.645 


± 





.058 


2 


.765 


± 





.057 





.792 


± 





.002 


6 


2881 


± 


20 


0.11 


± 


0.12 


0.15 


± 


0.17 


Gl 406 


7. 


.994 


± 





.015 


1 


.977 


± 





.015 





.794 


± 


0. 


.001 


6 


2887 


± 


20 


0.14 


± 


0.12 


0.18 


± 


0.17 


GJ 3147 


7 


.022 


± 


0. 


058 


2 


.773 


± 





.054 





.779 


± 





.002 


6 


2839 


± 


21 


0.16 


± 


0.12 


0.21 


± 


0.17 



Table 3 — Continued 



Name 




EW(N 


a 


I) 




EW(C 


a 


I) 




H 2 o-K2 


Sp.T. 


T £ 






[M/H] 


[Fe/H] 


LHS 292 


4. 


.791 


± 


0. 


.035 





.294 


± 





.035 





.759 


± 


0.002 


7 


2772 


± 


25 


-0.28 ± 0.12 


-0.41 ± 0.17 


Gl 644 C 


5. 


.117 


± 


0. 


.054 





.571 


± 





.051 





.745 


± 


0.002 


7 


2727 


± 


29 


-0.21 ± 0.12 


-0.32 ± 0.17 


GJ 1111 


5 


.661 


± 





039 


1 


.039 


± 





.038 





.758 


± 


0.002 


7 


2769 


± 


26 


-0.12 ± 0.12 


-0.19 ± 0.17 


LHS 2090 


6. 


.242 


± 





.032 


1 


.405 


± 





.032 


o 


.758 


± 


0.002 


7 


2769 


± 


26 


-0 03 + 12 


-0 06 + 17 


V* V492 Lyr 


6 


.406 


± 





.075 


1 


.225 


± 





.078 





.739 


± 


0.003 


7 


2709 


± 


28 


-0.02 ± 0.12 


-0.04 ± 0.17 


Teegardcns 


3 


.429 


± 


0. 


.039 





.269 


± 





.032 





.713 


± 


0.003 


8 


2637 


± 


30 


-0.38 ± 0.12 


-0.55 ± 0.17 


2MASS J1835 b 


3 


,886 


± 





.075 





.721 


± 





.073 





.692 


± 


0.003 


8 


2578 


± 


39 


-0.27 ± 0.12 


-0.40 ± 0.17 


LHS 2065 


7. 


.252 


± 


0. 


.104 


1 


.316 


± 





.099 





.687 


± 


0.003 


8 


2564 


± 


40 


0.12 ± 0.12 


0.16 ± 0.17 


LHS 2924 


4. 


,995 


± 





.049 





.845 


± 





.053 





.661 


± 


0.002 


9 


2492 


± 


43 


-0.13 ± 0.12 


-0.19 ± 0.17 


Gl 908 c 


2 


153 


± 





.115 


2 


.275 


± 





.113 


1 


.035 


± 


0.004 





3995 


± 


47 


-0.41 ± 0.12 


-0.59 ± 0.17 


Gl 809 c 


4. 


,085 


± 





.043 


4 


.193 


± 





.041 


1 


.057 


± 


0.002 











-0.15 ± 0.12 


-0.21 ± 0.17 


Gl 644 AB C 


3 


623 


± 





.048 


2 


.769 


± 





.045 


1 


.015 


± 


0.003 


1 


3845 


± 


19 


-0.27 ± 0.12 


-0.39 ± 0.17 


Gl 829 AB C 


4. 


,696 


± 


0. 


.123 


3 


.265 


± 





.129 





.907 


± 


0.004 


3 


3349 


± 


22 


-0.09 ± 0.12 


-0.13 ± 0.17 



a SPM J0011+5908 

b 2MASS J18353790+3259545 

c Stars with low quality K band spectra 



Note. - 
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Table 4. BT-Settl H 2 0-K2 indices 



T eff [M/H] 

-1.0 -0.5 0.0 +0.3 +0.5 



4400 


1.0633894 


1.0551468 


1.0459396 


1.0346350 






4300 


1.0610328 


1.0530623 


1.0447932 


1.0325842 






4200 


1.0555427 


1.0486027 


1.0440688 


1.0326834 






4100 


1.0496983 


1.0415743 


1.0399236 


1.0308537 






4000 


1.0419007 


1.0330761 


1.0336640 


1.0265948 


1.0123681 


3900 


1.0302266 


1.0223378 


1.0250625 


1.0202581 


1.0082819 


3800 


1.0083129 


1.0047255 


1.0108148 


1.0072626 


0. 


.99948109 


3700 


0.98369432 


0.98394169 


0.98987729 


0.99077235 





.98483467 


3600 


0.95944974 


0.95519382 


0.96753477 


0.96615374 


0. 


.96534533 


3500 


0.93851486 


0.93301470 


0.94241768 


0.93860547 





.94239536 


3400 


0.92003132 


0.91276339 


0.91955948 


0.91273634 





.91865987 


3300 


0.90412690 


0.89518129 


0.89743695 


0.89064011 


0. 


.89750592 


3200 


0.88679857 


0.87540384 


0.87727299 


0.86834456 


0. 


.88154377 


3100 


0.87051604 


0.85641549 


0.85691544 


0.84896337 


0. 


86984104 


3000 


0.85463719 


0.83051899 


0.83022602 


0.82209357 





.86153062 


2900 


0.84301951 


0.87462956 


0.80480058 


0.79553348 


0. 


.79291515 


2800 


0.83379220 


0.86284561 


0.77417570 


0.76562934 





.76116386 


2700 


0.83353582 


0.85048313 


0.74699251 


0.73209011 


0. 


.72973959 


2600 


0.83127154 


0.84092836 


0.71229056 


0.69373495 





.68974967 


2500 


0.84359432 


0.83413574 


0.68145572 


0.65795920 


0. 


.65185812 


2400 


0.83627202 


0.82175367 


0.64586819 


0.61983636 


0. 


.61785883 


2300 


0.84981127 


0.81318717 


0.62799732 


0.60425485 


0. 


.59720128 


2200 


0.85792452 


0.79609664 


0.58683221 


0.60165298 


0. 


.57924869 



Note. - 
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Table 5. M dwarf Metallicity Calibration Sample 



Name Sp.T. SPOCS Predicted 

[M/H] [Fe/H] [M/H] [Fe/H] 



Gl 212 


Ml 


0.16 


0.19 


0.02 


0.03 


HD 46375 B 


Ml 


0.20 


0.24 


0.21 


0.29 


Gl 872 B 


M2 


-0.16 


-0.22 


-0.17 


-0.25 


Gl 797 B 


M2 


-0.09 


-0.09 


-0.16 


-0.23 


Gl 250 B 


M2 


-0.01 


0.14 


0.01 


0.01 


NLTT 14186 


M2 


0.04 


0.05 


0.17 


0.24 


HD 222582 B 


M3 


-0.02 


-0.03 


0.02 


0.02 


Gl 324 B 


M4 


0.25 


0.31 


0.22 


0.31 


Gl 768.1 B 


M4 


0.12 


0.16 


0.08 


0.10 


Gl 231.1 B 


M4 


-0.08 


-0.04 


0.02 


0.02 


GJ 3348 B 


M4 


-0.10 


-0.22 


-0.01 


-0.02 


Gl 783.2 B 


M4 


-0.09 


-0.15 


-0.20 


-0.29 


Gl 544 B 


M4 


-0.15 


-0.18 


-0.05 


-0.09 


Gl 611 B 


M5 


-0.49 


-0.69 


-0.32 


-0.45 


Gl 777 B 


M5 


0.19 


0.21 


0.03 


0.03 


Gl 166 C 


M5 


-0.08 


-0.28 


-0.10 


-0.15 


NLTT 15867 


M5 


-0.05 


-0.10 


-0.14 


-0.21 


Gl 376 B 


M6 


0.11 


0.20 


0.11 


0.14 



Note. - 
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Table 6. Statistics of the [Fe/H] prediction models 



Calibration 


P 


N 


rms 


RMSp 


R 2 


Source 


B05 


5 


48 


0.483±0.047 


0.037±0.006 


0.858±0.028 


Table 


4« 


JA09 


1 


6 


0.082±0.016 


0.004±0.002 


0.472±1.224 


Table 


l 6 


SL10 


1 


19 


0.144±0.017 


0.022±0.005 


0.489±0.188 


Table 


l c 


Nil 


2 


23 


0.167±0.035 


0.034±0.010 


0.426±0.238 


Table 


2 d 



RA11 3 18 0.208±0.030 0.020±0.004 0.670±0.161 Tabled 



Note. — a) Bonfils et al. (2005), b) Johnson & Apps (2009), c) Schlaufman & 
Laughlin (2010), d) Neves et al. (2011) 
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Table 7. [Fe/H] Comparison Sample 



B05 JA09 SL10 W09 RA11 L07 
Name [Fe/H] [Fe/H] [Fe/H] [Fe/H] [Fe/H] Class 



Gl 699 






-0.76 




-0.39 


dM 


Gl 411 


-0.34 






-0.14 


-0.41 


dM 


Gl 905 




0.17 


0.14 




0.19 


dM 


Gl 447 






-0.28 




-0.01 


dM 


Gl 725 B 


-0.37 




-0.44 


-0.04 


-0.36 


dM 


Gl 725 A 


-0.32 




-0.29 


0.01 


-0.49 


dM 


GJ 1111 






0.25 




-0.19 


dM 


Gl 273 


-0.15 


0.07 


-0.06 


0.02 


-0.17 


dM 


Gl 860 AB 


-0.18 


0.03 


-0.10 


-0.15 


-0.11 


dM 


Gl 234 AB 


-0.00 


0.32 


0.22 




0.13 


dM 


Gl 628 


-0.12 


0.12 


-0.01 


-0.08 


-0.02 


dM 


Gl 687 


-0.13 


0.12 


-0.02 


0.04 


-0.09 


dM 


LHS 292 






0.20 




-0.41 


dM 


GJ 1245 AC 




0.20 


0.20 




-0.14 


dM 


GJ 1245 B 




-0.07 


-0.18 




-0.13 


dM 


Gl 876 


0.01 


0.38 


0.22 


0.02 


0.18 


dM 


Gl 412 A 


-0.52 




-0.50 


-0.33 


-0.40 


dM 


Gl 388 


0.03 


0.37 


0.20 


-0.21 


0.28 


dM 


Gl 873 


-0.10 


0.17 


0.03 


-0.21 


-0.01 


dM 


GJ 1116 AB 






0.89 




-0.12 


dM 


GJ 3379 


-0.12 


0.08 


-0.04 




0.08 


dM 


Gl 445 


-0.30 




-0.34 




-0.25 


dM 


Gl 526 


-0.26 


-0.05 


-0.17 


0.06 


-0.30 


dM 


LP 816-60 


-0.17 


-0.05 


-0.16 




0.06 




Gl 251 


-0.23 


-0.03 


-0.16 


-0.13 


-0.07 


dM 


Gl 169.1 A 


-0.01 


0.34 


0.20 


0.06 


0.36 


dM 


Gl 402 




-0.06 


-0.18 


-0.04 


0.20 


dM 


Gl 205 


-0.10 


0.12 




-0.13 


0.35 


dM 


Gl 213 




-0.12 


-0.23 




-0.25 


dM 


Gl 752 AB 


-0.03 


0.23 


0.09 


0.00 


-0.05 


dM 


Gl 285 


0.09 


0.55 


0.40 




0.40 


dM 
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Table 7— Continued 



B05 JA09 SL10 W09 RA11 L07 
Name [Fe/H] [Fe/H] [Fe/H] [Fe/H] [Fc/H] Class 



Gl 555 


0.00 


0.36 


0.23 


0.03 


0.22 


dM 


Gl 338 A 


-0.19 






-0.42 


-0.18 


sdM 


Gl 581 


-0.25 


-0.09 


-0.21 


-0.10 


-0.10 


dM 


Gl 338 B 


-0.34 






-0.29 


-0.15 


dM 


Gl 268 AB 


0.17 


0.76 


0.65 




0.10 


dM 


Gl 643 






-0.41 


-0.07 


-0.22 


dM 


Gl 1156 




0.07 


-0.01 




0.15 


dM 


Gl 625 


-0.48 




-0.50 


-0.35 


-0.23 


sdM 


Gl 408 


-0.22 


-0.02 


-0.14 


-0.09 


-0.09 


dM 


GJ 1286 




-0.01 


-0.10 




-0.04 


dM 


G 203-47 


-0.08 


0.19 


0.06 




0.33 


dM 


LHS 3376 






-0.68 




-0.12 


dM 


LHS 3799 




0.20 


0.14 




0.38 


dM 


GJ 1224 






-0.29 




-0.05 


dM 


LHS 3558 


-0.18 


0.04 


-0.10 


-0.14 


0.01 


dM 


Gl 686 


-0.39 




-0.33 


-0.25 


-0.28 


dM 


GJ 3146 










-0.11 


dM 


LHS 2065 










0.16 


dM 


Gl 849 


0.15 


0.53 


0.36 




0.31 


dM 


LHS 224 




-0.06 


-0.18 




-0.46 


dM 


LHS 1809 




-0.10 


-0.25 




-0.01 


dM 


LHS 3549 




0.05 


-0.05 




0.16 


dM 


Gl 176 


-0.06 


0.19 


0.06 


-0.16 


0.15 


dM 


LHS 495 




0.13 


0.03 




0.20 


dM 


LHS 252 




0.06 


-0.01 




-0.01 


dM 


V* V547 Cas 


-0.26 


-0.05 


-0.18 


-0.18 


-0.28 


dM 


Gl 436 


-0.03 


0.26 


0.10 


-0.05 


0.04 


dM 


GJ 1119 




0.14 


0.03 




-0.04 


dM 


Gl 649 


-0.15 


0.07 


-0.05 


-0.13 


-0.04 


dM 


GJ 3147 










0.21 


dM 


LHS 2924 










-0.19 


dM 
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Table 7— Continued 



B05 JA09 SL10 W09 RA11 L07 
Name [Fe/H] [Fe/H] [Fe/H] [Fc/H] [Fe/H] Class 



HIP 57050 


-0.02 


0.31 


0.15 


0.07 


0.05 


dM 


Gl 777 B 






-0.39 




0.03 


dM 


Gl 164 


-0.04 


0.23 


0.12 




0.02 


dM 


LHS 18 










-0.48 


dM 


Gl 179 


-0.01 


0.32 


0.16 


0.02 


0.23 


dM 


Gl 212 


-0.04 


0.18 


0.07 


-0.24 


0.03 


dM 


LHS 3577 


-0.26 


-0.05 


-0.18 


-0.22 


-0.11 


dM 


GJ 1214 




0.04 


-0.07 




0.20 


dM 


LHS 220 






-0.72 




-0.35 


dM 


LHS 3605 


-0.51 




-0.57 


-0.19 


-0.28 


dM 


Gl 611 B 










-0.45 


dM 


LHS 3593 




0.03 


-0.08 




0.03 


dM 


LHS 1706 




0.23 






0.21 


dM 


V* V492 Lyr 










-0.04 


dM 


GJ 3253 


0.05 


0.47 


0.38 




-0.07 


dM 


Gl 630.1 A 


-0.11 


0.15 


0.02 




-0.39 


dM 


HIP 79431 


0.19 


0.56 


0.39 




0.46 




V* V1513 Cyg 


-0.89 




-1.02 


-0.80 


-0.64 


sdM 


LHS 494 


0.08 


0.49 


0.32 




0.24 


dM 


NSV 13261 




0.31 


0.26 




-0.29 


dM 


G 168-24 


-0.29 




-0.26 


-0.09 


0.01 


dM 


LHS 3576 


-0.13 




-0.02 


-0.20 


-0.19 


dM 


LHS 3409 










-0.52 


sdM 


LHS 6007 


-0.04 


0.27 


0.14 




0.16 


dM 


G 210-45 


-0.16 


0.06 


-0.05 


-0.12 


-0.05 


dM 


HIP 12961 


-0.20 




-0.07 




0.01 




LHS 3591 


0.06 


0.35 


0.21 




-0.12 




G 262-29 


0.16 


0.38 


0.28 


-0.19 


-0.30 


dM 
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Table 8. M dwarf Wide Binaries [Fe/H] Comparison 







This Work 


B05 


JA09 


SL10 


W09 


Name 


Sp.Type 


[M/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 


GJ 1245 AC 


M6 


-0.09 


-0.14 




0.20 


0.20 




GJ 1245 B 


M6 


-0.09 


-0.13 




-0.07 


-0.18 




Diff: 




0.00 


0.01 




0.26 


0.3 




Gl 338 A 


MO 


-0.13 


-0.18 


-0.19 






-0.42 


(jl 338 B 


MO 


-0.11 


-0.15 


-0.34 






-0.29 


Diff: 




0.02 


0.03 


0.15 






0.13 


Gl 725 A 


Ml 


-0.34 


-0.49 


-0.32 




-0.29 


0.01 


Gl 725 B 


M3 


-0.25 


-0.36 


-0.37 




-0.44 


-0.04 


Diff: 




0.09 


0.13 


0.05 




0.15 


0.05 


V* V547 Cas 


M2 


-0.19 


-0.28 


-0.26 


-0.05 


-0.18 


-0.18 


LHS 115 


M3 


-0.19 


-0.27 










Diff: 




0.00 


0.01 










Gl 412 A 


Ml 


-0.28 


-0.40 


-0.52 




-0.50 


-0.33 


Gl 412 B 


M6 


-0.27 


-0.39 










Diff: 




0.01 


0.01 










Gl 643 


M3 


-0.15 


-0.22 






-0.41 


-0.07 


Gl 644 C 


M7 


-0.21 


-0.32 










Diff: 




0.06 


0.10 
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Table 9. Best M dwarf Targets for Planet Searches 



Name 


PMSU 


Sp.T. 


V 


K s 


[M/H] 


[Fe/H] 


H a 


v sin i 2 


Pop 


Gl 581 


2420 


M2 


10.57 


5.84 


-0.06 


-0.10 


-0.359 


<2.1 a 


Thin 


Gl 687 


2797 


M3 


9.15 


4.55 


-0.06 


-0.09 


-0.34 


<2.8 a 


Thin 


Gl 408 


1709 


M2 


10.03 


5.50 


-0.06 


-0.09 


-0.307 


<2.3 a 


Thin 


Gl 251 


1094 


M2 


9.9 


5.3 


-0.04 


-0.07 


-0.272 


<2.4 a 


Thin 


Gl 752 AB 


3037 


Ml 


9.12 


4.67 


-0.03 


-0.05 


-0.42 


<2.5 b 


Thin/Thick 


Gl 649 


2673 


Ml 


9.7 


5.6 


-0.02 


-0.04 


-0.49 


<1.9 C 


Thin 


Gl 628 


2599 


M3 


10.1 


5.08 


-0.01 


-0.02 


-0.234 


l.l d 


Thin 


Gl 447 


1849 


M4 


11.12 


5.65 


-0.00 


-0.01 


-0.181 


<2 C 


Thin 


LHS 1809 


1004 


M4 


14.5 


8.4 


-0.00 


-0.01 


-0.231 


4.3 f 


Thin 


LHS 252 


1355 


M6 


15.05 


8.67 


0.00 


-0.01 


0.312 


4 f 


Thin/Thick 


LHS 3558 


3203 


M3 


10.4 


5.9 


0.01 


0.01 


-0.039 


<3.2 a 


Thin 


G1250 B 


1092 


M2 


10.1 


5.7 


0.01 


0.01 


-0.292 


<2.5 b 


Thin 


Gl 212 


956 


Ml 


9.9 


5.8 


0.02 


0.03 


-0.391 


<2.5 b 


Thin 


Gl 436 


1830 


M3 


10.7 


6.1 


0.03 


0.04 


-0.37 


<l c 


Thin/Thick 


LHS 3549 


3195 


M5 


14.0 


8.1 


0.12 


0.16 


-0.131 


<4.5 f 


Thin/Thick 


Gl 905 


3743 


M4 


12.3 


5.9 


0.14 


0.19 


0.092 


<1.2 a 


Thick 


Gl 876 


3604 


M2 


10.2 


5.0 


0.14 


0.19 


-0.2 


<2 C 


Thin 


Gl 402 


1687 


M3 


11.7 


6.4 


0.15 


0.20 


-0.28 


<2.3 C 


Thin 


Gl 555 


2310 


M3 


11.3 


5.9 


0.16 


0.22 


-0.177 


2.7 c 


Thin 


Gl 324 B 


1383 


M4 


13.2 


7.7 


0.22 


0.31 


-0.249 


<2.5 b 


Thin 


Gl 849 


3478 


M2 


10.4 


5.6 


0.23 


0.31 


-0.384 


<2.5 b 


Thin 


Gl 205 


930 


MO 


8.0 


4.0 


0.25 


0.35 


-0.401 


l d 


Thin/Thick 


Gl 169.1 A 


780 


M3 


11.1 


5.7 


0.26 


0.36 


-0.174 


1.9 C 


Thin/Thick 





2 Reference of H Q 



Gizis et al 



(2002) 



1992 


), (d) 


Reinert 


2007 


), (e) 



Delfosse et al 



( 1998bJ), (b) 



Mohantv fc Basril |2003J), and (f) 



Browning et al 



(2010 ), (c) 



Jenkins et al 



(2009) 



Marcv fc Chen 



Table 10. M dwarf Planet Hosts and Metallicity Calibration Sample 





This Work 




KHM 


B05 


JA09 


W09 


SL10 


R10 


M sin(i) 




Name 


Sp. Type 


[M/H] 


[Fe/H] 


Sp. Type 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 


Mj 


Planet Notes 


HIP 79431 


Ml 


+0.33 


+0.46 


M3 


+0.16 


+0.52 




+0.35 


+0.60 


2.1 


Jupiter 


Gl 849 


M2 


+0.23 


+0.31 


M3 


+0.14 


+0.58 


+0.20 


+0.41 


+0.49 


0.82 


Jupiter 


Gl 179 


M3 


+0.17 


+0.23 


M3.5 




+0.30 


+0.02 


+0.20 




0.82 


Jupiter 


GJ 1214 


M4 


+0.15 


+0.20 


M4.5 




+0.03 




+0.28 


+0.39 


0.0179 


super-Earth 


Gl 876 


M3 


+0.14 


+0.19 


M4 


+0.03 


+0.37 


+0.02 


+0.23 


+0.43 


2.64 


2 Jupiters + Neptune + super-Earth 


Gl 176 


M2 


+0.11 


+0.15 


M2.5 




+0.18 


-0.16 


+0.06 




0.0265 


super-Earth 


Gl 436 


M3 


+0.04 


+0.05 


M2 


-0.03 


+0.25 


-0.05 


+0.10 


-0.00 


0.072 


Neptune 


HIP 57050 


M3 


+0.03 


+0.04 


M4 


-0.02 


+0.32 


+0.07 




+0.12 


0.298 


Neptune 


HIP 12961 


Ml 


+0.01 


+0.01 


MO 








-0.07 




0.35 


Neptune 


Gl 649 


Ml 


-0.02 


-0.04 


Ml 


-0.18 


+0.04 


-0.13 


-0.03 


+0.14 


0.328 


Neptune 


Gl 581 


M2 


-0.06 


-0.10 


M3 


-0.25 


-0.10 


-0.10 


-0.22 


-0.02 


0.0492 


Neptune + 3 super-Earths 



